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Abstract 


A model for the prediction of wake related sound generation 
by a single airfoil is presented. The basis of the model is 
an assumption that the net force fluctuation on an airfoil may 
be expressed in terms of the net momentum fluctuation in the 
near wake of the airfoil. This near wake model results in the 
forcing function for sound generation being the spectra of the 
two point velocity correlations in the turbulent region near 
the airfoil trailing edge. 

The spectra of the two point velocity correlations were 
measured for the longitudinal and transverse components of 
turbulence in the wake of a 91,4 cm (36 inch) chord NACA series 
63-009 airfoil. The airfoil was installed in the 7 by 10 tunnel 
at the NASA Ames Research Center. The test velocity was 45.7 
m/sec (150 ft/sec) and incidence angles were set at 0° and 6°. 

Two x-probe hot wire sensors were used in the measurements. 

Cross power spectra of the signals from these sensors provided 
the data required for the two point velocity correlations. Some 
preliminary experiments were conducted using a 25.4 cm (10 inch) 
chord airfoil similar to the larger one in a free jet facility 
at the University of Washington. 

A scaling procedure was developed using the turbulent bound- 
ary layer thickness as the primary measure. The model was then 
used to predict the radiated sound from a 5.1 cm chord airfoil 
for which acoustical data were available. Agreement between the 
predicted and measured sound radiation spectra was good both in 
level and spectral shape. 

The single airfoil result was extended to a rotor geometry 
and comparisons studied for various aerodynamic parameters. 


The motion of the source was found to influence the spectxmn 
of radiated sound very little. This was due to the relatively 
broad band spectrum of the source. Sample computations were 
made for large subsonic rotors and the noise levels due to the 
wake mechanism alone was found to be significant. This was 
due in part to the concentration of the acoustic energy in the 
frequency range 2000 to 9000 Hz which is the range of highest 
subjective response by the human ear. 
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1 . Q Introduction 

The noise radiated from helicopter rotors and other air- 
craft fans consists of a coinbination of pure tones and broad 
band effects. The pure tones typically occur at frequencies 
which are the blade passage of the rotor and its higher 
harmonics. When the Mach number of operation approaches unity 
other pure tones, associated with the rotational frequency 
appear. These are the "multiple pure tones" or Buzz Saw com- 
ponents" which one hears when tbe current super jet takes off. 
The blade passage frequency pure tone mechanism appears to be 
due to the steady state pressure distribution on the individual 
blades and the interactions between the rotor and stator rows. 
Recently, there has been evidence to suggest that the interac- 
tion of inflow turbulence with the rotor blading will give rise 
to this pure tone as well. 

Broad band noise is usually associated with the inter- 
action of the blading with turbulence in the flow, either with 
inflow turbulence or wake turbulence of the airfoils themselves 
Possible mechanisms of this type include boundary layer turbu- 
lence, random static pressure fluctuations related to the wake 
shedding at the trailing edge of the blade, radiation from the 
stresses in the wake turbulence itself; and force fluctuations 
caused by the interaction of the airfoil with inflow turbulence 
All these mechanisms are present and there are probably more 
that are yet undefined. The one which dominates the others de- 
pends upon the flow being considered and the geometry of the 
body in the flow. There is ample evidence to suggest that the 
dominant mechanism is dependent upon both the Reynolds and Mach 
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numbers of the flow. 

The present study is an attempt to relate noise generation 
by a single airfoil operating in a smooth turbulence free flow 
field to the properties of the turbulence in the wake of the 
airfoil. Then, once the mechanism is identified, to extend the 
model to a rotor to provide a design tool. What is hoped from 
this approach is that a model for the prediction of the spec- 
trum of broad band noise generation in the design phase of a 
turbomachine may be developed. 

Extension of the results from the single airfoil to a 
rotating array of airfoils was made analytically using an exist- 
ing multipole analysis. By using the single airfoil results 
as a source of information on the statistical properties of the 
fluctuating forces, the effect of rotation of an array of blades 
could be determined. 

It is expected that for large rotors such as may be expect- 
ed to be used in large STCL and helicopter aircraft, the 
blades and hence their wakes will be large enough that wake 
related noise will be an important feature in the noise consid- 
erations of these designs. Because the wake is related to 
aerodynamic design and certain possibilities exist for its 
modification (for example, boundary layer suction to reduce the 
size of the turbulence correlation volumes) it is important to 
understand the detailed mechanism so that flow control can be 
directed most efficiently. In current designs, the small size 
and high speed of operation make the frequency of wake related 
noise too high to be of practical interest (except for some 
tip vortex shedding on large helicopter rotors). This also 
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tends to degrade the usefulness of typical model rigs for study- 
ing this mechanism; they are too small in size. However, when 
Mach numbers are reduced to decrease the tone noise and blad- 
ing size is increased to accomodate the aerodynamic loading, 
then the wake noise will be in a frequency range of very criti- 
cal subjective importance. 

The most important wake parameters in terms of noise gen- 
eration appear to be the actual size of the wake and the velo- 
city of operation. Normally one sees the sound power expressed 
as proportional to velocity to the 5,5 or 6th power. The effect 
of wake correlation volume is as important as velocity since 
the acoustic power scales with the geometric size (correlation 
volume times physical volume) scale to the sixth power. As a 
result, one might experience a sixth power variation nf acoustic 
power with geometric scale. Unfortunately it isn't that simple 
and there appear to be few shortcuts to the rather pragmatic 
procedure of comparing designs one by one using the most com- 
plete model available. It is only in this way that all the 
effects of a given change can be properly accounted for. Gen- 
eral dependencies can of course be stated but when the actual 
practical designs are put down it will be a complex interrela- 
tionship of parameters which will determine the correct choice. 
For an example, a designer might feel that by reducing the velo- 
city of operation he can reduce the sound generation; but, when 
he designs the fan he finds he must increase the chord to main- 
tain the same aerodynamic loading. When he tests the machine 
he finds it to be as noisy but with the sound shifted to another 
part of the spectrum. In order to design for noise one must 
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have a clear understanding of the implications of design changes. 

In this study an attempt is made to relate sound, genera- 
tion to wake properties so that the effects of such things as 
geometric scaling can be assessed. Experimental measurements 
of the cross power spectra of two point velocity correlations 
were made since these were found to provide the forcing function 
information needed by the model. The results were extended to 
a multib laded rotor so that design considerations could be 
made using a model more closely approximating the actual rotor. 
An effort has been made to simplify the final results to a 
point where they may be used by the designer who may not be 
a specialist in the details of aerodynamic noise theories and 
methods . 


2.0 Review of the Literature 


This review will be restricted to the literature which 
applies most directly to the problem of sound generation by 
bodies immersed in smooth turbulence free flows. This means 
the concern will be with those investigations which have attemp- 
ted to relate sound generation to either the static pressure fluc- 
tuations on the surfaces of the body or to the turbulence which 
is generated by the airfoil as the flow passes over it. Pub- 
lished efforts dealing with airfoils are somewhat limited and 
for the most part recent since motivation for this research has 
been public concern with aircraft noise. Some earlier work has 
relevance. This is the study of Aeolian tones generated by 
flows over cylinders. The very first works were concerned with 
musical instruments and later with loadings on telegraph lines 
and other cables. This early literature will be referred to 
only in the context of its relevance to the problem of aero- 
dynamic noise and how it relates to aerodynamic performance 
parameters. 

Modem studies of aerodynamic sound generation are based 
on solutions of the non-homogeneous wave equation developed by 
M. J, Lighthill (1,2,3). It is his work which has provided the 
basis for the large volume of analysis in the field of aero- 
dynamic sound generation. By combining the equations of con- 
servation of linear momentum with that of conservation of mass, 
Lighthill was able to write a non-homogeneous wave equation 
which would govern the radiation of aerodynamic sound from a 
region of turbulent motion. Justification of the approach is 
best expressed in Lighthill' s own description. It is important 
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to understand this reasoning to best judge its scope of appli- 
cability. To quote Lighthill: 

Considering a fluctuating fluid flow occupying a limited 
part of a very large volume of fluid of which the remain- 
der is at rest. Then the equations governing the fluc- 
tuations of density in the real fluid will be compared 
with those which would be appropriate to a uniform 
acoustic medium at rest, which coincided^ with the real 
fluid outside the region of flow. The difference be- 
tween the two sets of equations will be considered as 
if it were the effect of a fluctuating external force 
field, known if the flow is known, acting on the said 
uniform acoustic medium at rest, hence radiating sound 
in it according to the ordinary laws of acoustics. 

In this way all the forcing effects are lumped into a forcing 
function from which the resulting sound radiation may be pre- 
dicted. The consequence of the method is that one must be 
able to specify the unsteady aerodynamics of the source region. 

To show the way the model is constructed, consider the 
equation of conservation of linear momentum and use the summa- 
tion convention where repeated indices imply summation; 


at 



and the equation of conservation of mass: 





Taking the partial derivative of continuity with respect to 
time and the partial derivative of linear momentum with respect 
to the space variable then subtracting the two equations, one 
can express the result as a wave equation of the form 
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where T. . is the “Lighthill tensor" and is defined as 
ij 

Tij =^UiUj - oij - 0.1 pSi^ . 


In these expressions ; pO ~ density, a^ = speed of sound in the 
fluid at rest, * = stress tensor, u. = velocity in direction 

1. J 3- 

x^ (i = 1,2,3). It is this wave equation which forms the basis 
for the many theories recently developed for sound generation 
by aerodynamic flows. 

The general Kirchhoff solution can be used to solve the 
non-homogeneous wave equation (Sobolev, reference 4). This 
original theory by Lighthill did not include the effect of 
solid bodies contained in the region of fluid motion. This was 
done by N. Curie (5). Curie wrote the Kirchhoff solution for 
the case of aerodynamic flows as 





r 


4 


t 

4tt 



or 


dr 


* SJ l ^ 1 ^ 

an ^ ^aoT dn at 



where- r =j‘x - *y | . 

The integrands of both terms are taken at the retarded time. 
This simply means in the present case that the sound observed 
in the far field was generated at the source location at an 
earlier time and the relevant source strength is that which 
occurred at the earlier time. The first integral is over the 
volume of the turbulent region and the second is over the sur- 
face of the active region. 

In studying the solution for the case of immersed solid 
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bodies, Curie was able through, several manipulations to cast 
the solution in a particularly revealing form'. He was able to 
show an equivalent solution as 



Here again the solution is the sum of a volume and a surface 
integral. The sound produced is due to a volume distribution 
of quadrapoles represented by the fluctuating stresses and a 
surface distribution of dipoles represented by the pressure on 
the surface of the body. is exactly the force per unit area 
exerted by the solid boundaries on the fluid in the x^ direction. 
This second Integral is the additional effect of the solid 
boundary to the original solution written by Lighthill. The 
solution presented by Lighthill was for no boundaries and is 
exactly the first integral of Curie's result. 

Using dimensional analysis, Curie developed a relationship 
between sound power and aerodynamic parameters. The total acous- 
tic power output was shown to take the form 

where U is a typical velocity and L is a characteristic length. 

Another approach to the solution of the wave equation with 
immersed bodies was presented by P. E. Doak using the method of 
Green's functions (6). Once the appropriate Green's function 
is known for a space, the solution can be Written immediately. 
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Actuall}'', the Green’s function approach has not been used much 
because free space problems have proven difficult enough. In 
worh on tone noise generation the method has been used. Doak 
wrote expressions for the radiated sound intensity in terms of 
the covariances of the first time derivative of the surface 
pressure. It was these forms which provided the starting point 
for I. J, Sharland's (7) prediction of noise generation by 
small flat plates immersed in jet flow and the extension to 
broad band sound generation by fans. For a flat surface with 
dimensions small with respect to the wavelength of the radiated 
sound, the total radiated sound energy could be written: 



where t' is the retarded time, the surface S is in the Xj^, X2 

(or XpX^) plane, p is the pressure difference fluctuation at 

(xj^,X2) and the overbar denotes a time average. The integration 

is taken over only one side of the plate. Since the wavelength 

is larger than the plate dimension, xj^ and x^ can be taken as 

equal to Xj^ and X2* To express the integral, a proportionality 

factor S can be used which has the units of area and is called 
c 

the correlation area. The integral may then be written as 



S 


and the total radiated acoustic energy as 
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Sharland now related the pressure to a fluctuating lift force 
per unit area. He wrote the pressure derivative as 

where U is the local mean velocity parallel to the plate and 
Cj^ is described as a lift coefficient. 

Sharland used some boundary layer assumptions to express 
the vortex shedding noise and derived the expression for over- 
all sound power; 


E = -?-10 CU dK 
0 ?„ 1 


Span 

He compared this prediction with some measurements of overall 
sound pressure for one inch chord fiat plates located in the 
potential core of a one inch diameter round jet. His radiation 
pattern suggested a dipole as did the power law of the sound 
power with velocity. A dipole yields a dependence upon velocity 
to the sixth power. A plot of the radiated sound as a function 
of velocity taken from Sharland* s paper is shown in Figure 1. 

It may be seen that the agreement is quite good. This is extra- 
ordinary in view of the simplicity of both the model and the 
experimental measurements. The agreement does provide encourage- 
ment that the general approach is valid and further work along 
this path would be worthwhile. 

The efforts of Sharland were the first to attempt the pro- 
blem of wake related sound generation from airfoils, but there 
were some some attempts using circular cylinders. The most 
complete approach using modem theories was made by 0. M. 

Phillips on vortex sound generation from small cylinders (8). 







Radiated SPL at 90° to Plate, (dB re to 20 N/m 


LI 



Velocity, 


ft/sec 



Figure 1 Sound Radiated from a One Inch Flat Plate 
Immersed in a Smooth Flow 
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Fluctuations in lift and drag were computed for Reynolds num- 
bers from 40 to 160 (based on diameter) using the velocity field 
measurements of Kovasznay (9). Kovasznay bad measured the vel- 
ocity field downstream in the near wake of a circular cylinder 
and found he could express it in the functional form 



Phillips wrote the drag and lift forces using the momentum eqxia- 
tion (integral form) for incompressible flow. For the drag 



Drag could be computed directly from the velocity field data 
but the lift computation would require the normal component of 
turbulence which had not been measured. Phillips used the vor- 
ticity equation along with continuity to obtain the lift fluc- 
tuation. He obtained 


L — D COS Simt 


The axial length scale was found to be very large for Re 100 
and for the range 100*<Re<160 to be about 17D (D = cylinder 
diameter). Using this^ the radiated dipole sound was expressed 


~ o.2Tcos®0 


a§r^ 


where 1 = cylinder length, D = diameter, and S == Strouhal num- 
ber (a/ 0.2). 


When higher Re 5 niolds numbers were tested a similar expres- 
sion was derived but vji.th. a numerical factor of 0.037 rather 
than the 0.27 used in the above. Also, the span correlation 
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was found to decrease to a value of about 3D at the highest 
Reynolds numbers tested. 

Measurements which gave extensive data on the Strouhal 
number have been presented both by Kovasznay and by G. H, 
Gerrard (10). Garrard measured lift and drag fluctuations on 
circular cylinders using a pressure transducer mounted on the 
top on the cylinder. He concluded that the drag fluctuations 
are an order of magnitude smaller than those of the lift. The 
radiated sound intensity is given without derivation as 


X ^3 ^2 U O 

2 0 5 

where (C^) * is the lift coefficient which was determined ex~ 
perimentally to be 0. 8 at Re = 7 x 10^ and 0.01 at Re = 4 x 10^. 
It appears that the expression is incorrect as written since 
there is no provision for total cylinder length. Probably the 
1 should be replaced by Sl^, the product of the correlation 

U. Q 

length and the cylinder length. The trend with Reynolds number 
of the lift coefficient seems reversed to this writer since one 
would expect a higher intensity at the lower Reynolds number. 
This is due to the decrease in the correlation length with an 
increase in Reynolds number. 

The chronology of the above research was such that each 
investigation was done at a different time so that each had 
access to the results of the previous investigators. The de- 
velopment is perhaps more systematic than that which has been 


done more recently when the various workers were proceeding in 
parallel . 

The next three programs to be discussed were ail directed 
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to the airfoil wake noise problem. Each of the three was done 
in a different part of the world and was done independently of 
the others. It is interesting to compare the approaches as 
they are completely different although they started at a common 
point which was the techno logy level in 1965. 

I. Kavrak used a hydrogen bubble flow visualization tech- 
nique (11) in water to observe the wake of both circular cylin- 
ders and airfoils. Through still photographs and films, he 
was able to qualitatively estimate the flow field. Based upon 
these observations, he asstomed pressure waveforms and computed 
the spectrum of the static pressure in the wake. Kavrak noted 
that vortices shed from the airfoil were neither regular in 
size or shape. They tended to grow while being convected down- 
stream but did not change shape. The vortex cross section ap- 
peared to be elliptical with an axis ratio ranging from 2 to 6. 
The vortex pattern for an airfoil wake was like that of a rect- 
ified signal, not thealternating vortex pattern typically seen 
in flow around cylinders. 

In a later paper (12), Kavrak measured the fluctuating 
static pressure in the wake using a probe designed for the pur- 
pose. He observed the most intense pressure fluctuations in 
the regions of maximxam shear stress. The observations were used 
to develop a model based on curve fits to the various measured 
parameters. The vortices were assumed to be elliptic waveshapes. 
Pressure-time histories were fitted using least squares methods 
and the resulting functions were Fourier transformed to yield 
spectra of the pressure. Once the various- curve fits required 
were completed, the maximiam pressure in the wake was expressed 
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in tile functional form 



v/here F(Re) is a Reynolds number dependent function, 0( ~ angle 
of incidence, = incidence angle at minimum drag, b Is the 
transverse distance from the center of the wake to the 1/2 max- 
imum velocity point. Sound measurements were made with a rotor 
made up of four airfoils located with their spans parallel to 
the axis of rotation. A microphone placed on this axis measur- 
ed the radiated sound , Comparisons of wave shapes between the 
radiated sound and the pressure in the wake were made. No act- 
ual quantitative comparisons were presented comparing source 
pressure to radiated sound. levels. For such a geometry, one 
would expect significant interaction effects due to the wakes 
acting as an inflow turbulence mechanism. The frequency content 
of these interactions may be such as to avoid the problem. 

While the above research was being carried out in Istanbul, 
B. D. Mugridge was approaching the problem somewhat differently 
in England. Mugridge (13) attempted to relate sound generation 
to the covariance of static pressure on. one side of an airfoil. 
He obtained a relationship for the power spectrirn of the force 
fluctuation on the blade and identified the fluctuations with 
eddy shedding. Assuming the eddies to be convected along at 
some fraction of the free stream velocity, he was able to ex- 
press the force spectrxmi for the airfoil. In any method of this 
type, some measure of the spatial extent of the correlation must 
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b© made. In this case it is the correlation area on the air- 
foil surface. The method used by Mugridge was to compare the 
radiated sound to the auto spectrum of the static pressure at 
one point on the airfoil surface. By forming a ratio of these 
two quantities in a way which would account for the geometry of 
the radiation path, an estimate of the correlation area could 
be obtained. The idea was to measure the cause (the surface 
pressure) and the effect (the radiated sound), then determine 
the correlation area required to get the correct answer. The 
experimental measurements were made on a custom airfoil design 
(which appears to have a thickness to chord ratio of 0,20) 
mounted in a wind tunnel. The airfoil was cambered so that the 
boundary layer turbulence would be most confined to the suction 
surface. This would validate the one sided integration proce- 
dure. Radiated sound was taken as that picked up by a pressure 
transducer mounted on the wall of the wind tunnel. Area ratios 
based upon the ratio were 

A^/A = 1/5 at 1 kHz, and = 1/45 at 2.5 kHz 

When compared with results taken by Schloemer (14) for boundary 
layers with adverse pressure gradients, some estimate of what 
might expected was obtained. These estimates were 

A /a = 1/15 at 1 kHz, and A /A = 1/90 at 2.5 kHz 
No explanation was advanced for the differences. Mugridge con- 
cluded the problem would be resolved when experimental measure- 
ments of A could be made directly using two surface mounted 
static pressure transducers. The measurements were made with 
a cambered airfoil and efforts were made to insure the bulk of 
the turbulence would be generated on the suction side of the 
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airfoil. One way to eliminate the integration on the back side 
is to make the boundary layer very thin on that side. That was 
apparently the aim in Mugridges’s method. 

In a later paper (15), Mugridge presented cross correlation 
measurements between the surface static pressure at 84% of the 
chord with the longitudinal component of turbulence in the boun- 
dary layer near the trailing edge of the airfoil. He found sub- 
stantial correlation between the surface static pressure and 
the velocity in the outer edge of the boundary layer. This 
correlation had not been seen in similar measurements taken on 
wall boundary layers. The large eddy structure in the outer 
region was felt to be relatable to the higher pressure levels 
measured for airfoil surfaces as compared to those pressure 
levels measured in wall boundary layers at similar flow condi- 
tions. The airfoil surface pressure levels were typically 15 dB 
above those for the wall boundary layer. Obviously there was 
something significant about the trailing edge. 

Some measurements of overall turbulence length scales were 
made in the wake of the airfoil using spatial correlations. 

These were found to collapse fairly well if the wake momentiom 
thickness was used in the normalization. 

The third approach to the airfoil sound generation problem 
was made by the present author (16), In this program the net 
force fluctuation on the airfoil was related to the momentum 
fluctuations in the airfoil wake. The relation was established 
using the unsteady form of the integral momentum equations. 

The component of momentum normal to the airfoil chord was taken 
as the forcing function for the sound generation. This required 
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measurement of the transverse component of turbulence. In this 
model the correlation required was the volxime of correlated 
momentum in the wake of the airfoil. This was obtained through 
measurements of the correlation lengths in the three coordinate 
directions. The correlation voliime was taken as that defined 
by the two point velocity correlations measured using two sin- 
gle wire sensors. This meant that the volume for the transverse 
component would be assumed to be identical to that for the 
stream component. The spectrum of the normal component was 
determined using an x-probe hot wire sensor. In this way even 
if the volume was determined using the longitudinal component 
the magnitudes would be correctly measured. These results were 
then combined to predict the spectrum of the radiated sound. 

The resulting expression was 



When this prediction was ‘compared with measured sound rad- 
iation data for the same airfoils, fair agreement was obtained. 
The experimental sound radiation data was obtained for the air- 
foils located in the potential core of a free jet discharging 
through an anechoic chamber. In the facility both the magnitude 
and the radiation pattern could be obtained. The predicted 
levels were usually lower than those measured over the entire 
frequency range. The prediction has now been seen to contain 
an error in the formulation. The same mistake was made as that 
mentioned before by Gerrard. The dependence on the span corre- 
lation should be on the product of the span correlation length 


J. 
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and tinie actual span length not on the square of the span corre- 
lation length as was shown above. Fob the particular model 
above and the data with which it was compared, the difference 
would add 6 dB to the entire spectral level predicted. This 
increases the predicted spectra to levels having better agree- 
ment with those measured. 

The major failing of the model was the poor prediction of 
the spectral levels at higher frequencies. There was a consis- 
tent tendency to overpfedict these levels. This was known to 
be due to the usage of a constant value for the correlation 
lengths where in fact the correlation length is expected to de- 
crease with increasing frequency. It was the development of 
this more complete model where the correlation lengths would be 
measured for the particular velocity component of interest and 
the lengths would be determined as functions of frequency which 
formed the basis for the present study. 

Measurements of radiated sound from single airfoils have 
been recently presented by a group at the United Aircraft 
Research Laboratories (17), This laboratory has built a very 
complete facility for the measurement of sound radiation from 
airfoils. The facility consists of a large open throat wind 
tunnel having its test section in an anechoic chamber. At the 
time of the present report the background noise levels were 
quite high, hopefully these will be reduced to levels allowing 
better resolution in test results. In the data reported, a 
pure tone was observed over a wide range of Reynolds numbers, 
in fact several pure tones were observed at various operating 
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conditions. No explanation other than vortex shedding was ad- 
vanced. This tone phenomenon was observed to be dependent upon 
the incidence angle, increasing the incidence angle increased 
the probability of appearance of the tone(s). Surface pressure 
correlations in the fom of two point space correlations were 
made using transducers. In the fully turbulent, no tone, re- 
gion no correlation was observed in the span direction. The 
transducers used were 1/2 inch diameter microphones located on 
a 9 inch chord airfoil. The diameter of the transducers is 
larger than the correlation length. Because of this the auto- 
spectra taken must be held suspect since the response is surely 
an average of the pressure over the face of the microphone. 

This might be expected to substantially underestimate the sur- 
face static pressure. 

It was because of the small correlation lengths which this 
author had observed in previous work (16) that the approach of 
relating the net force fluctuation on the airfoil to the momen- 
tum fluctuation in the wake was taken. The greater spatial re- 
solution and ease of probe placement versus the problems asso- 
ciated with positioning and calibration of surface transducers 
was felt to justify the more indirect wake momentum approach. 
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3.0 Models for the Prediction of Sound Radiation 

In this section the models used for prediction of sound 
radiation by both single airfoils and rotors are outlined. The 
motivation for the development of a single airfoil model is 
twofold. As a component of a rotor, understanding the sound 
generation characteristics of the airfoil can help in understan- 
ding similar mechanisms on the rotor. Secondly, if sufficient 
confidence can be established in the concept of relating net 
force fluctuation to near wake turbulence, then this approach 
may be useful in studying noise generation by the airplane it- 
self exclusive of the engine noise. 

The rotor model is based upon a rotating array of fluctua- 
ting forces simulating the forces on the blades . The force on 
each blade is defined in terms of the experimental measurements 
and the single airfoil model. 

3.1 Single Airfoil Model 

It may be seen from the Curie solution to the wave equation 
that three possible mechanisms exist for the generation of sound 
by a stationary airfoil in a smooth turbulence free flow. These 
are the boundary layer related pressure fluctuations, force 
fluctuations due to some net effect like vortex shedding, and 
turbulent stress fluctuations of the boundary layer and wake 
turbulence. The second mechanism has been chosen as the domi- 
nant one for the present case. Boundary layer noise has been 
discounted based on results of both Sharland and Mugridge which 
indicate boundary layer noise to yield too small a sound level 
(down approximately 15 dB from measured results. Direct radia- 
tion from wake turbulence has been discounted because of 
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measured radiation patterns and velocity dependencies of SPL, 
Experimental results for several surface contours, incidence 
angles, and Mach number range indicate a dependence of overall 
radiated SPL on velocity to the 5.5 to 6th power. The examples of 
these dependencies may be seen in Figures 2 and 3. Direct rad- 
iation from the turbulent stresses would result in an 8th power 
dependence on velocity. The fact that the sound radiation is 
observed to occur at the same frequency as the transverse com- 
ponent of turbulence in the wake further supports the dominance 
of this dipole mechanism. This similarity in the spectra may 
be seen in the results presented in reference 16. 

If the sound generation mechanism is accepted as some form 
of surface wake interaction, what can be done to relate it to 
aerodynamic performance? Recall the far field solution for a 
net force fluctuation on a body immersed in a flow was 

The most natural approach might be to form the density covar- 
iance then observe that the forcing function is the covariance 
of the surface static pressure. The primary reason that this 
approach was not taken in the present investigation was the con- 
cern that the correlation lengths required for determination of 
the correlation areas would be smaller than the spatial resolu- 
tion capabilities of the transducers which would be available 
for the measurements. It was thought that if a method relating 
sound generation to the aerodynamics of the wake directly could 
be established then perhaps a simpler relationship for the mec- 
hanism in terms of normal aerodynamic performance would result. 
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In addition, the measurement of the turbulence structure of the 
wake could be carried out using hot wire sensors thereby allow- 
ing reasonably good spatial resolution of the correlation 
lengths . 

Since the sound generation has been shown to be relatable 
to the net force fluctuations on the airfoil, a simple method 
can be used to relate these force fluctuations to wake turbu- 
lence. If the Gauss divergence theorem is applied to the above 
integral for the pressure on the body and the momentum equation 
is used to relate the resulting pressure gradient to velocity, 
the following correspondence can be written:* 



This states that the net force fluctuation is equal to the net 
momentum fluctuation in the airfoil wake. 

The idea of the model is the wake immediately following 
the airfoil will reflect the net momentum on the blade since it 
is here that the boundary layers from both sides will mix, sta- 
tic pressure variations on both sides will balance and the re- 
sultant effect will be reflected in the momentum in the wake in 
the region of mixing. If one takes surface integrals on one 
side of an airfoil and then predicts the sound generation based 
upon that, in effect, it is assumed that the two sides are un- 
correlated. It seems to be more reasonable to allow the effects 
of the two sides to average out in the near wake then use that 
averaged effect as the measure of the net force fluctuation on 


^Details are presented in Appendix A. 
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on the airfoil. 

The net force flnctnation on the airfoil is then to be 
approximated by the momentiira in the wake volxmie immediately 
downstream of the airfoil. Details of the development are some- 
what lengthy and are presented in detail in the appendix. The 
near wake is taken as the correlation volume nearest the trail- 
ing edge of the airfoil. This volume is defined as the product 
of the correlation lengths in the three coordinate directions. 
The integration for the resulting covariance of the two point 
velocity correlations would be the wake thickness, the airfoil 
span, and the total correlation length in the stream direction. 
The particular correlation was found to be somewhat different 
for the longitudinal and transverse components of turbulence. 
This is expected since the integral scales for the longitudinal 
and transverse components are different even for isotropic tur- 
bulence; see for example the discussion in reference 24. 

Using the above assumptions, the spectrum of sound radia- 
tion can be related to the spectra of two point velocity corre- 
lations as 

2 2 2 

ya Ys 

With: 

47ra^/x/ ^ 4^7ra^/xl 
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and"^ is the angle of the observer relative to a normal to the 
chord of the airfoil, Figure 4. 

The variables L^, L 2 , and are correlation lengths de- 
fined as the frequency dependent integral scales in the three 
coordinate directions. These lengths are obtained from spatial 
integrations of a coherency function for the two point velocity 
correlations. In the experimental part of the program measure- 
ments of L^(uj^u^,6o) and L^(u2U,2j<a 3) were made directly but be- 
cause of the amount of data reduction the lengths for the cross 
terms were not done. These lengths were taken as the arithme- 
tic mean of the lengths determined for the longitudinal and 
transverse components. They were approximated by 

Li(UiU^,Oi) = L^(u 2 U^, to) = 1/2 (L^(u^u^,& 3 ) + L^(u 2 U 2 ,cO )) 

i = 1,2,3 

There is contribution from the cross terms for non-zero inci- 
dence angles. When the wake is not symmetric, the integrals 
containing the cross terms do not vanish. For the data obtain- 
ed in this Investigation the cross terms made very little con- 
tribution to the radiated sound. 

There results four terms in the expression for the radiated 
sound: 

(l)i=l, j=l corresponds to a dipole oriented along 
the stream direction. This term has been neglected 
in previous analyses based upon observations that 
the alternate shedding of vortices from opposite sides 
of the body would give a quadrapole in the stream 
direction. This appears to be valid for low Reynolds 
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number flows where pure vortex shedding exists ; how- 
ever, for high Reynolds number flows and particularly 
with non-zero incidence angles it appears that longi- 
tudinal momentum fluctuations will be present. 

(2) i=2, j==2 corresponds to the fluctuation in the 
transverse direction due to the transverse momentum. 
This is the mechanism used in previous models of vor- 
tex sound generation. It will yield a dipole aligned 
in the direction of a normal to the airfoil chord. 

(3) i=2, j =lori=l, j=2 are the cross terms which 
are present when the airfoil operates at non-zero in- 
cidence angles. Because of the sign on the radiation 
angle these two terms will always subtract from the 
sum of the other two terms. The effect is to render 
the radiation pattern non- symmetric. 

Contributions of these various effects will be illustrated 
in the results section of this report. In this section, the 
way the model includes the effects of non-zero incidence angles 
will be shown along with other aerodynamic effects. This basic 
model forms a starting point for the development of an engineer- 
ing model which will eventually allow scaling effects to be de- 
termined. The extension of these results to the rotor is the 
next subject to be considered here, 

3.2 Rotor Model 

The analytical model chosen for extension of the single 
airfoil results to a rotor geometry was that presented by 
Ffowcs Williams and Hawkings (21). They solved the non-homo- 
geneous wave equation for an array of sources (simple, dipole, 
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and quadrapole) moving in a circular path. Details of the de- 
velopment of their analysis into a model for the present appli- 
cation are presented in Appendix B. 

The equation to be solved is 

It has the general solution 





which may also be written in terms of its generalized Fourier 
transform: 



For the case of point forces, the source term reduces to 




X 






k.^ ajxl 



where is the magnitude of the dipole force. 

These general results are reduced for the particular case 
of an evenly distributed array of blades each with its force 
uncorrelated with the others. The forces on the blades are 
assumed to be stationary random processes. With these assump- 
tions, the radiation of sound from a multibladed rotor can be 
approximated as 
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[fj ,4-nxi)d^ (e^, nii-i} 

n^-oo te=i 

V J 

where d^ is the component of the force in the direction of the 
observer. 

Notice that the assumption that the entire blade span may 
be approximated by a single force vector is quite restrictive 
in that the blades of practical interest are usually long and 
the aerodynamic properties vary significantly along the span. 

A span distributed result can be obtained by breaking the rotor 
into annular segments then adding the sound radiation which re- 
sults from each segment to get the total radiation from the 
rotor. The reason this simple procedure worked is because the 
correlation which exists in the span direction is small. Each 
radius along the blade span radiates essentially independently 
of the others so the use of a simple summation is valid. 

3.2.1 Relation to Single Airfoil Results 

The next step is to use the single airfoil wake model to 
describe the force on the individual blades of the rotor. The 
total force on the blade is taken as the integral of the pres- 
sure over the blade surface. This is the same procedure used 
in the single airfoil model development. The force is 

d(6j,,-b)=r -pty,^;) dscv*) 

From the single airfoil model write 

j yty.t) ds{y) = f ^ p'^Ly,-^)^y 

5 ■'y 






y.. 
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Now transform to express this in terms of frequency: 






The result for the modulus of the force fluctuation is 
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The results have been left in terms of the power spectrum of 
the two point velocity correlations. These correlations can 
now be expressed within the assumptions of the single airfoil 
model as 

/ f 
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This result was programmed to form a prediction of the radiated 


sound from a rotor due to the random force fluctuations exper- 
ienced during operation in a smooth turbulence free flow. 


3.2.2 Computer Program for Rotor Sound 

Discussion of the computer program developed for the rotor 
sound prediction will be brief; the approach was fairly straight 
forward. The complexities in coding the model results in a 
fairly substantial effort to get the computer program checked 
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out. 

Characteristics of the sound mechanisms are known relative 
to the individual blades as given by the single airfoil model. 
These blades are arrayed to form a rotor; each blade is set at 
an angle to the rotor axis to properly align with the relative 
velocity vector (staggtr angle plus incidence). In order to 
express the properties on the individual blades in the rotor 
coordinate systenij two changes of basis were required. The 
first was to unrotate the stagger angle and the second was to 
go from the blade to the observer coordinate system. 

Details of the algebra are not presented here. It is suf- 
ficient to state that the program was checked out using constant 
force vectors on the blades and comparing the resulting tone 
noise generation with other available methods. Another check 
of the method was made by setting the rotor speed to zero so 
that the rotor degenerated to a set of stationary airfoils. 

The spectra for this case agreed with the single airfoil model. 
It was felt the rotor program was functioning correctly. 
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4.0 Experimental Program 

The experimental program consisted of two parts. Initially 
tests were made using a 25.4 cm chord airfoil located in the 
potential core of a free jet tunnel at the University of 
Washington. Correlations of wake turbulence using two single 
sensors were made to provide some information on the physical 
dimensions of the correlations to be measured in the larger 
test facility at the NASA Ames Research Center. These measure- 
ments in the Ames 7 by 10 subsonic wind tunnel #1 using a 91.4 
cm chord airfoil will form the major part of the experiemntal 
results . 

The purpose of the experiment conducted at Ames was to de- 
termine the power spectra of the two point velocity correlations 
in the wake of a single airfoil. The airfoil used was repre- 
sentative of contours, used in helicopter rotors (NACA 63-009). 

A large airfoil, 91,4 cm chord, was used to insure the measure- 
ments would be made at a high Reynolds number and thereby avoid 
the uncertainties existant in the transitional flow regime. 

Estimates of the two point velocity correlations, u^u^, 
^1^2’ ^^2^1’ ^2^2 made using two x-probe hot wire sen- 

sors. Cross power spectral densities of signals from these 
four hot wire sensors were used to obtain the required spectral 
data. Notice that correlations of the velocities in the span 
direction were not studied. The two point correlations in the 
span direction would be expected be strictly a volume turbulence 
effect since there is no surface for them to exert a force on. 
They are parallel to the airfoil. A single traverse was made 
in order to estimate the magnitude of this turbulence component 
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and as expected it was found to be comparable to the other two 
components. 

As a check of the accuracy of the four wire method addi- 
tional cross power spectra were taken using two single wire 
sensors. This yielded the spectrum. It was felt that 

data from these single wires are probably more accurate than 
those from the x-probes because of the reduced possibility for 
probe interference effects. It was expected this increased 
scope in the experiment would help in evaluating the confidence 
limits of the four wire method. 

Data reduction was performed from data stored on magnetic 
tape. Details of the test facility, instrumentation, and test 
procedure are given below. 

4.1 Test Facility 

As mentioned above the Ames 7 by 10 subsonic xi?ind tunnel 
#1 was used for this investigation. This tunnel is normally 
used to measure forces and moments needed to evaluate aerodyna- 
mic performance of flow models. The tunnel is a recirculating 
type driven by electric fans. Speeds in excess of 60 m/sec are 
readily available, although in this test a fixed speed of 45.6 
tn/sec was used. The entire two story building in which the 
test section is located is sealed x^hile the tunnel is running 
which means the control room pressure is equal to the static 
pressure in the test section. Leaks do not cause any particu- 
lar difficulty with this arrangement since there is little or 
no flow through them. Control of the test section velocity is 
simply done using a rheostat. The tunnel held its setting quite 
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well for most of these tests. However, when the airfoil was 
set to the 6 . degree angle of attack, some drift was experienced. 

The airfoil used was a 91.4 cm (36 inch) chord NACA series 
63-009 contour having zero camber angle. The aspect ratio of 
this airfoil was approximately 1.5 so that care had to be exer- 
cised to insure the flow was straight and parallel. Tufts were 
attached along both sides of the airfoil above and below the 
projected traverse plane and the flow was observed to be accept- 
ably parallel in the region of the measurements when the angle 
of incidence was at zero degrees. Only one survey point was 
taken at 6° incidence because it was expected that the flow 
would be three dimensional and unreliable for an airfoil of 
such a low aspect ratio . 

Velocity in the test section was set using a pitot static 
probe mounted in the test section upstream of the model. Read- 
out was made on a vertical water manometer. Pressure variation 
was less than 2.5 mm of water for the tests reported here. 

Turbulence levels of the test section were measured with 
the movable probe traversed out into the potential flow. The 
level was found to be very low. The longitudinal component 
was measured to be 0.5% on an RMS basis and the transverse com- 
ponent was less than this. The measurements were made by re- 
cording the peaks of the autocovariance functions and the levels 
were too low for the transverse component to load the SAICOR 
Correlator used for the measurement. The level could have been 
determined using a true RMS meter had it been convenient and 
necessary to the purpose of the experiment. 

The Reynolds number corresponding to the test airfoil and 
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flow condition was approximately 2 x 10°(91,4 cm chord and 
45.6 m/sec velocity). 

Two x-type hot wire sensors were used, one stationary and 
the other movable relative to the stationary one. The station- 
ary probe was attached to a floor mounted strut which had an 
elliptical cross section. The photograph in Plate I shows the 
installation. The stationary probe used a right angle adapter 
in order to get the correct wire orientation while the movable 
probe was on a sting mount allowing it to be positioned very 
close to the stationary probe. The movable probe was mounted 
on the tunnel traverse mechanism. This mounting allowed all 
three coordinate directions to be traversed by remote control. 
Backlash problems required that probes be always positioned from 
one direction. This was done throughout the testing and offered 
no particular operational difficulties. Estimated position un- 
certainties for the three directions were: 
stream; -.01 cm + .01 cm 
span: -.012 cm +.012 cm 

normal: -.012 cm +.012 cm 

Backlash in the span and normal directions was about .05 cm but 
this error was avoided by the test procedure described above. 

The stream direction traversing was done with a DISA traversing 
mechanism which had no discernible backlash. 

4.2 Instrumentation 

Data flow for the experiment was from the hot wire ane- 
mometers to magnetic tape then later from magnetic tape to the 
time series analyzer. A list of the instrumentation is given 
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below and a flow chart of the data flow is shown in Figure 5. 

Hot Wire Anemometry and probes: 

constant temperature anemometer DISA 55D01 

linearizer unit DISA 55D10 

stationary x-probe DISA 55A39 

movable x-probe DISA 55A38 

single wire probe DISA 55A22 

DC voltmeter DISA 55D30 

Tape Recorder: Ampex FR-1300A, 14 channel 

IRIG single extended, bandwidth 10 kHz at 30 ips 

Time Series Analyzers : 

Correlator: SAICOR model SAI-43A 

Time Series Analyzer: Time Data model 1923 

4.3 Test Procedure 

Test procedure consisted of calibrating the probes, set- 
ting the wire sensitivities, positioning the probes and record- 
ing the data on magnetic tape. More details follow. 

4.3.1 Probe Calibration 

Calibration was carried out in the test section by locat- 
ing the probe in the free stream away from the airfoil and vary 
ing the velocity over a range of test speeds anticipated (24 
m/sec to 48 m/sec). The exponent required to set into the lin- 
earizer was then determined and set into the linearizer to 
yield an output proportional to velocity in the test range. 

The effectiveness of the linearization was checked by perform- 
ing another calibration on the linearized output. The response 
were found to be linear over the test velocity range. 

4.3.2 Operation 

Typically the sensitivity of a hot wire will vary a small 
amount from wire to wire. In this test the gain on the linear- 
izers was adjusted so that all the wires were equal sensitivity 



FIGURE 5 SCHEMATIC FOR DATA ACQUISITION AND REDUCTION 
OP TWO POINT VELOCITY CORRELATIONS 
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Initially the wires were set at the same overheat ratio and 
temperature compensation probes were used on all four hot wire 
circuits so the adjustment to achieve equal sensitivity was 
very small. This adjustment was done by locating the movable 
probe in the free flow region and setting the linearizer gain 
so that the DC levels of each of the two wires was equal. Then 
this probe was moved into position directly downstream of the 
stationary probe, that is, so the x-arrays were located in 
line at an axial separation of approximately 2 mm. The DC lev- 
els of the stationary probe were then set to match those of the 
movable one. In this way all four wires tvould respond to the 
same steady velocity with the same voltage. We observed that 
moving the downstream probe spanwise until it was clear of the 
wake of the other probe made no difference in the DC levels 
and concluded that the procedure was sufficient for setting 
probe sensitivities equal. When making correlation volume mea- 
surements, the movable probe was placed manually close to the 
fixed probe. The procedure used was to position the movable 
probe directly behind the fixed probe which had been located 
on the center line of the wake and a given distance do^TOstream 
of the trailing edge. Spacing between the centers of the 
x-probes was set at approximately 2.5 mm. From this position, 
traverses were made in the three coordinate directions. Actual 
selection of spacing was made through the on-line use of the 
SAICOR correlator. With the cross correlation displayed on the 
scope, we could select separations which would provide a smooth 
decrease in the correlation between the two probes. If this 
on-line capability were missing, positioning would have been 
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very much more trial and error. 

Once the probes were positioned, a sample of data was re- 
corded on magnetic tape. Sample lengths varied from one to two 
minutes. Data recorded were: 

®A1’ ^A2’ ®A1 ®A2’ ®A1 " ®A2' 

®B1" ®B2’ ®B1 ®B2’ ®B1 " ®B2 

where , 

= voltage from probe A wire 1 
e ^2 voltage from probe A wire 2 
= voltage from probe B wire 1 
6^2 = voltage from probe B vjire 2 
These data could then be analyzed to yield the cross power 
spectral densities of the two point velocity correlations. 

Data were recorded at four axial positions for the station- 
ary probe and an array of correlations were made at each of 
these positions. The axial positions used were: 3.5, 8.3, 16.7, 
and 36% of the chord downstream of the airfoil trailing edge 
and on the center line of the wake. The wake centerline was 
located as the point of minimum mean velocity. 

The two probes were positioned as shown in Plate II. As 
mentioned earlier, the steady flow direction was perpendicular 
to the stationary probe body. This is not an ideal condition, 
but for these correlation measurements there was little choice 
if close spacings between the two probes were to be obtained 
without wake interference problems. One difficulty with such 
a mounting method is the strain gauge effect in the wire resis- 
tance which results as an effect of vortex shedding induced 
forces on the wire supports. If this occurs the supports 
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vibrate and one receives false signals at frequencies near 
15 kHz. For the present case these signals were outside the 
frequency range of interest and masked by the levels of wake 
turbulence. Correlations taken on data where the probe was 
outside the wake did indicate the vortex shedding mentioned. 

The tape recorded was band limited at 10 kHz so this did not 
appear in the recorded data. 

It was observed that one very good way for locating the 
wake centerline would be to use the zero sheaf stress since 
this was available on-line as the cross correlation of the 
"sum" with the "difference" signals taken on a single x-probe. 
4,3.3 Data Reduction 

Data reduction was carried out on two different time series 

analysis instrimients . The on-line instrument, the SAICOR, was 

also used for data reduction of the final recorded data. This 

was primarily for the single x-probe data to determine the 

wake property distributions and cross correlations of the data 

from the two single wires to be used as a data check of the 

four wire method. The settings used on the SAI-43A were: 

Full linear averaging with 128 summations of 1024 points 
each. The work position was varied in order to fill the 
frame on the X-Y plotter. Very smooth correlations 
resulted from this instrument used in this way. 

The cross power spectral densities were performed directly on 

the Time Data 1923 using its Fast Fourier transform capability. 

By cross correlating the sum and difference signals from the 

two probes, the spectra of the two point velocity correlations 

were obtained. The form used was the magnitude of the power 

spectrum. 
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The forms may be summarized as: 
u^(y,o) u^*(y,£0) 
u^(y,tt>) U 2 *(y,«) 

U2(y,«) u^*(y,w) 

U2(y,w) U2*(y,uJ) 

An alternate procedure to the usage of the Time Data Analyzer 
was to directly Fourier transform the correlation functions ob- 
tained on the SAI-43A correlator. This was done using a time 
series analysis program developed as part of this program. The 
program developed would give cross power spectra, autospectra, 
and present magnitude and phase information. 

For data presentation the various plots were digitized 
using a graphics machine at the tfW Computer Center which allows 
plots to be accurately transcribed to binary tape. The trans- 
cribed data was then conditioned to the proper physical units 
using further programs and plotted using a CalComp Plotter. 

A parabolic fit package was used for the correlation functions 
and a simple straight line plot was used for the spectra. 

4,4 Sources of Errors in Measurements 

Errors in the experimental data are larger than had been 
expected. The most significant errors are not due to the un- 
certainties in probe placement as had been expected or in the 
accuracy of the reading, but to what are thought to be operator 
errors in the data reduction and probe interference errors in- 
herent in the two probe correlation method. 

As was mentioned earlier, data reduction was carried out 
using both a Time Data Analyzer and a SAICOR Correlator. Re- 
duction using the Time Data Analyzer was not performed by the 
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author but by a professional operator employed for that purpose. 
This was the only way that data could be reduced using the ma- 
chine. Reduction on the SAIGOR was done by the author. We 
were never able to make the Time Data results agree quantita- 
tively with either the SAICOR results or with meter readings 
of the true RMS levels. The RMS levels computed from the 
correlation functions reduced on the SAICOR agreed very closely 
with those measured using true RMS meters. 

Reduction was carried out using the Time Data Analyzer on 
two different occasions and significantly different results 
were obtained. When comparisons were made between spectra re- 
duced by Fourier transforming the covariance functions obtained 
on the SAICOR analyzer and those done on the Time Data Analyzer, 
a disagreement of 10 dB was observed. This was for the first 
set of data reduced. When the second set was done disagreement 
of 3 dB between the two sets of data, both reduced by the same 
operator on the same machine, occurred. 

As a result of the uncertainty in the actual levels, an 
alteimate procedure was used to determine the correlation length 
estimates from this data. The cross spectra were normalized 
with the autospectra of the velocity at the stationary probe. 
This formed a kind of coherency function. Spatial integration 
of this function yielded the correlation lengths as a function 
of frequency which was the proposed use of the data. All data 
finally used in the airfoil wake model were reduced in the first 
set of Time Data results. These are thought to be self consis- 
tant. The same gain error appears to be present in all this 

The normalized data should not reflect the gain setting 


data. 
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error. The results are in good agreement with the correlations 
measured using the SAICOR derived spectra. 

Another uncertainty in the data appears to be due to inter- 
ference between the probes. Notice in Figure 25 that the coher- 
ence function falls off rapidly as the probes are very close 
to each other then as the spacing increases the correlation in- 
creases again. This effect appears to be due to some sort of 
probe interference. It seems to be present for the stream cor- 
relations where the probes are located downstream of each other 
and the effect was not observed when measurements with two sin- 
gle wire sensors were used to repeat the correlations. The in- 
tegrations tended to reduce the error in terras of final effect 
on the correlation length. A similar interference probably 
exists for the transverse component but there was no way to 
check the x-probe measurements for this component. Future mea- 
surements of this type might be better done using two single 
wires slanted at 45° to provide the various components rather 
than the x-probe s used here. The time for the experiemnt would 
be four times as long but the results should be sufficiently 
more accurate to justify the method. The best estimate of the 
error in the correlation lengths at each frequency is 15% based 
on comparisons between the single wire measurements and the 
x-probe results. 

In the results presented here, the Time Data values are 
corrected for the suspected 10 dB error by normalization with 
the autospectra. For the SAICOR reduced data, in that excellent 
agreement was obtained with various data checks it is thought 
that these data are correct. 
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Finally, it should be noted that these measurements are 
somewhat involved and time consuming and errors of the magni- 
tudes noted appear to be somewhat common in this type of data. 
For the purposes of the acoustic predictions, the expected 
maximum scatter is i 3 dB which is sufficient although not as 
good as would be desired. 

4.5 Experiments with a 25.4 cm Airfoil 

Some preliminary testing was done using a 25.4 cm chord 
airfoil at the University of Washington. This airfoil was also 
an NACA series 63-009 contour. The test facility used was a 
free jet of rectangular cross section 20.4 cm by 45.7 cm. Sha- 
ping of the nozzle was done so that the wall curvature decreased 
uniformly to zero curvature at the exit plane thereby insuring 
straight parallel flow in the jet discharge. The airfoil was 
mounted in the potential core of the jet. 

Correlation measurements were made using two single wire 
sensors. Space time correlations for the longitudinal compon- 
ent of turbulence was measured in the three coordinate direc- 
tions and at four positions in the wake. 

General comments regarding technique and data reduction 
are similar to those for the 91.4 cm airfoil except that all 
reduction was done using a SAICOR model 43A correlator. The 
primary purpose of this experiment was to establish a test plan 
for the program using the large airfoil. As it turned out 
the data were able to be used in the scaling law development 
described later in this report. 
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5.0 Experimental Results 

The primary experimental results obtained in the program 
are the measurements of the wake turbulence for the 91.4 cm 
chord NACA series 63-009 airfoil. Test conditions common to the 
data presented for this airfoil are sumarized below: 

Airfoil type: NACA series 63-009, 91,4 cm (36 inch) 

Total temperature: 65° F 

Static pressure: 1 atmos. 

Velocity, free stream: 46.2 m/sec (150 ft/sec) 

Incidence angles: 0° and 6° 

5.1 Overall Wake Properties 

Wake surveys were made at five downstream positions using 
the movable x-probe hot wire sensor. Signals from this sensor 
were reduced to yield mean velocity, the longitudinal component 
and the transverse component of turbulence. 

Mean velocity profiles are presented in Figure 6, RMS tur- 
bulence profiles for the longitudinal component in Figure 7, 
and for the transverse component in Figure 8. These distribu- 
tions were used in the airfoil prediction model- In addition, 
some data were checked with longitudinal turbulence component 
measurements made with a single wire sensor. RMS levels between 
the two probes, the x-probe and the single wire, agreed to 
within 6%. 

A mean velocity profile was taken 1.2 mm downstream of the 
airfoil trailing edge. This was obtained using a single wire 
probe so that very fine spatial resolution would be obtained. 

The profile is shown in Figure 9. This profile was used to com- 
pute the wake momentum and displacement thicknesses used in the 
data evaluation. Drag coefficients based on the measured momen- 
tum thickness agreed with published results very closely. This 
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helps provide confidence in the measured profile. 

An indication of the wake development can be seen from the 
centerline changes of the properties along the stream direction. 
Centerline distributions of the mean velocity and the two turbu- 
lence components are shown in Figures 10 and 11. Notice that 
beyond 0.10 chord downstream the mean velocity defect is de- 
creasing at a 0.5 power as would be expected for a wake in the 
self preservation range. In terms of displacement thickness 
this corresponds to 70 thicknesses. Reference to turbulence 
literature suggests that self preservation does exist by this 
distance, but only for the mean velocity not higher moments. 

5.2 Power Spectra of the Turbulence Components 

The power spectra for the longitudinal and transverse tur- 
bulence components were obtained by a Fourier transform of co- 
variance functions determined using the SAICOR correlator. 

These power spectra are presented for each of the five stream 
positions for the wake centerline. This centerline spectrum 
was observed to be representative of the turbulence in the main 
part of the wake. 

Power spectra for the longitudinal component are presented 
in Figures 12 through 16. Notice that there appear be three 
spectra on each figure. The data shown are a Fourier ti"’ .sform 
of both sides of the measured auto covariance function and a 
cross spectrum of the entire function (two sided). This method 
gives a measure of the error in the entire procedure used in 
the data reduction from the auto covariance measurement through 
its manual digitization and Fourier transform. Notice that 
errors are negligible except for the high frequencies where the 
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effects of the noise floor are becoming apparent. 

The bandwidths (uncorrelated Fourier coefficients) for the 
data are 520 Hz. This tends to flatten the spectra in the low 
frequency range. The Time Data results have a bandwidth of 41 
Hz so better describe the details of the low frequency spectrum. 

Notice in the spectra for the longitudinal component that 
there is mono tonic decrease in the level with frequency. For 
the transverse turbulence, Figures 17 through 21, the energy 
does not exhibit this behavior but peaks at a higher frequency 
then decreases beyond this point. This is the general behavior 
noted by Tennekes and Liunley (24) for isotropic turbulence. 

Autospectra for the two turbulence components for the 6° 
incidence angle are presented in Figures 22 and 23. The major 
effect of the increase in incidence angle is to shift the ener- 
gy to lower frequencies as would be expected. 

5.3 Spectral Distribution of the Two Point Velocity Correlations 

The purpose of the measurements of the spectra of the two 
point velocity correlations was to determine the spatial corre- 
lation of the various frequency components. 

A measure of the spatial correlation was required as a 
function of frequency for the two turbulence components in each 
of the three coordinate directions. The way this measure was 
established was to normalize the magnitude of the cross power 
spectra with the autospectra for the signal from the stationary 
probe. This resulted in a spatial correlation function at each 
frequency and in each coordinate direction. This was essential- 
ly a cross plot of the space time power spectra. A spatial in- 
tegration of each of these frequency dependent correlation 
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Figure 23 Auto-Spectrum for the Transverse Component 
on the Wake Centerline, i = 6° . v,/C = .03 
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functions would yield an integral length scale as a function of 
frequency . 

The cross spectra data will be presented in the form of the 
spatial correlation functions at each of the frequencies. This 
cross plotting procedure is somewhat lengthy so only enough 
data will be shown to illustrate the method. 

The spatial correlations for the longitudinal component at 
0° incidence are shown for the stream direction in Figure 24, 
for the span direction in Figure 25, and for the normal direc- 
tion in Figure 26, Correlation plots for the 6° incidence case 
are presented in Figures 27, 28, and 29 in the same order as 
the 0° incidence data. 

The spatial correlation functions for the transverse com- 
ponent are sho^-m in the three directions in Figures 30, 31, and 
32 for the 0° incidence data and in Figures 33, 34, and 35 for 
the 6° case. 

The numerical integration of plots of this type was car- 
ried out on the CDC 6400 to yield a spectral distribution of 
the integral lengths in the different directions. Because of 
the few points in the spatial coordinates, the method used for 
the integration was the trapezoidal rule. A more sophisticated 
quadrature formula appeared to introduce more uncertainties be- 
cause of the scarcity of data paints. 

5.4 Frequency Dependent Correlation Lengths 

The correlation lengths measured through the above describ- 
ed integration for the stream direction are shown in Figures 36 
and 37. Data are presented for both components of turbulence 
and each of the two incidence angles studied. Notice that the 
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Component - Incidence =6° 
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Figure 30 Coherence Function in the Stream Direction - Transverse 
Component - Incidence =0° 
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Figure 33 Coherence Function in the Stream Direction - Transverse 
Component - Incidence =6° 
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lengths for the 6° case are larger at lower frequencies but 
smaller at the higher frequencies for the longitudinal compon- 
ent. For the transverse component the lengths are about the 
same up to 1 kHz where the 0“ lengths become as much as twice 
those for the 6° data. 

In the span direction the correlation lengths of both 
components are nearly the same over the entire frequency range 
for the 0° data, Figures 40 and 41. The same general behavior 
noted for the longitudinal component can also be seen here. 

The lengths are about the same at low frequency for both inci- 
dence angles with the 6^"* data rolling off at a lower frequency 
than the 0° data. 

The greatest differences occurred with the normal direction 
data. Substantially larger correlation lengths were observed 
in the 6^ data for both components of turbulence. Ratios of 3 
to 4 times the 0“ lengths were observed at different frequencies 
as may be seen in Figures 38 and 39. Notice that the normal cor- 
relation length for the 0° case was found to be .15 cm for the 
transverse component; this seems very small for a wake of this 
size. 

In general, the lengths measured for the transverse com- 
ponent are substantially different in the three coordinate di- 
rections. The extent of correlation in the stream direction is 
more that ten times that in the normal direction. An eddy length 
of 2 cm does not seem physically unreasonable for this flow 
since one half the wake width is about this size. The 0.15 cm 
normal length seems small; however, examination of both the cross 
power spectra and the cross correlation functions support this 


length as correctly measured. 

A rough check of the general magnitude of the stream cor- 
relation length can be made through a comparison with some eddy 
size estimates. Those suggested by Tennekes and Lumley (24) 
will be used. On page 259 of this reference a measure of the 
eddy size is taken as 

Ig = 2-ir/k 

where k is the wave number which is defined for the present 

data as k=cO/u with U the convection speed in the centerline 
o o 

region of the wake. One half this eddy size, to be consistent 
with the definition of the correlation lengths, is plotted in 
Figure 38. As may be seen the eddy size estimate is of the same 
order as the correlation lengths plotted but with substantial 
spectral differences. 

Another measure which can be made is to compute the inte- 
gral scale for each of the directions using the spatial corre- 
lation functions which are the values of the correlation func- 
tion at zero time lag. These results are shown in Figures 38, 
39, 40, and 42 since these were the only cases for which spa- 
tial correlation functions were available. The spatial integral 
scale is close to the zero frequency value of the correlation 
lengths. This is expected since the zero wave number for a 
generalized Fourier transform does correspond to the integral 
scale. 

The greatest differences occurred with normal direction 
correlation lengths. Substantially larger lengths were obser- 
ved in the 6° data for both components of turbulence. This 
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Figure 36 Frequency Dependence of the Correlation Length 
in the Stream Direction - Stream Component 
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increase in the correlation length across the wake might be ex- 
pected since the wake itself is thicker. The boundary layer 
on the suction side of the airfoil is thickened due to the in- 
creased loading. An increase of four times seems a large amount 
but not unreasonable since at 6“ the lift coefficient is approx- 
imately 0,6 which represents a resonably high loading on this 
airfoil. In the current data this was observed at the rather 
high Reynolds number of 2 x 10 so it's doubtful it is a transi- 
tion effect. 

Correlation lengths measured for the longitudinal component 
are shown in Figure 42 for the 25,4 cm chord airfoil at 0° inci- 
dence angle. These will be used in the scaling model. 

5.5 Reynolds Number Dependence of Correlation Lengths 

Three different airfoils which had basically the same con- 
tours were used to provide some infomation on the Re 3 molds num- 
ber dependence of the correlation lengths. The correlation 
length, based upon the flat part of the spectrum, to chord ratio 
for the three airfoils is shown as a function of Reynolds number 
in Figure 43, The smallest airfoil was a 5.1 cm chord series 

5 

65-010 operated at Re = 2 x 10 , the second was a 25.4 cm chord 

5 

series 63-009 operated at 3 x 10 , and the third was the 91.4 cm 

fi 

chord series 63-009 operated at 2 x 10 . Distributions of the 
correlation lengths for the longitudinal component in the flat 
part of the spectrum are plotted for the three directions. No- 
tice the normal and stream correlation length ratios decrease 
reasonably slowly with increased Re)niolds number in a way which 
might be expected since the -.2 power line plotted in the fig- 
ure is the same as expected for flat plate boundary layer growth 
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dependence on Reynolds number. The major difference occurs 
with the span correlation length which decreases much more ra- 
pidly in the transistional region than in the fully turbulent 
region. There appears to be a very strong dependence of the 
span correlation on the Reynolds number for airfoils of this 
t)^e. In fact, based on this plot, there could be a ten fold 

increase in the span correlation length over the range of 

5 6 

Reynolds number from 2 x 10 to 2 x 10 . Certainly much more 
work is needed to establish this dependence. That it is impor- 
tant is evidenced by the fact that if all else were equal the 
change in correlation length observed would result in an in- 
crease of 10 dB in the radiated sound pressure level. 

5.6 Phase Measurements in the Near Wake 

Phase angle measurements for different axial separations 
were determined as a function of frequency. This was done by 
computing the magnitude and phase representation of the cross 
power spectra. Results for the near wake are shown in Figure 
44, The phase angle may be seen to plot as a straight line for 
each of the separations. Computations between the frequency 
and separation indicate the same convection velocity for all 
the frequency components. This suggests that the turbulence 
can be considered as being converted along at the same mean 
velocity in a kind of frozen structure. 

Similar plots were not made for the transverse component 
but examination of the cross covariance functions suggested 
that the same general result would be valid. It appears that 
there may be some slight irregularities in that there is a 
small amount of distortion of the correlation function as it 
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convects downstream. This would suggest different velocities 
for different frequency components. The effect is not large, 
5.7 Normalized Autospectrum 

In an effort to compare some of the current data with that 
of other investigators a normalized plot of the transverse com- 
ponent of turbulence is shown in Figure 45. The free stream 
velocity and the wake displacement thickness was used to nor- 
malize the frequency. The resulting plot seems to put the 
energy in the same frequency range as Mugridge found for the 
static pressure fluctuations. 




6.0 Discussion 


Tlie results of the turbulence measurements were sufficient 
to complete the prediction models for the single airfoil and 
the rotor. In this section, the results of the predictions for 
the single airfoil will be compared with some experimental data 
in an effort to assess the effectiveness of the model, A scal- 
ing procedure is outlined and some sample computations for a 
large subsonic rotor design are presented. 

6.1 Single Airfoil Results 

The first question to be addressed is certainly how well 
does the model work. There are very few acoustic data avail- 
able for airfoils operating in smooth flows where wake effects 
may be assumed to dominate the sound generation. The best that 
could be done was to scale the present data to the range of 
parameters used in the acoustic tests reported in reference 16. 

A summary of the pertinent variables follow: 

Turbulence Measurements Sound Measurements 

series 63-009 airfoil series 65-010 airfoil 

45.7 m/sec 88.2 ra/sec 

91,4 cm chord 5.1 cm chord 

Re = 2 X 10® X 10^ 

incidence = 0° incidence - 0'^ 

As may be seen the Reynolds numbers differ by an order of 

magnitude. The basic contours are similar; the 65 and 63 series 

airfoils are not substantially different. Overall RMS wake 

turbulence spectra were available for the 5.1 cm airfoil and 

these were used for this prediction. Correlations were taken 

from the 91.4 cm airfoil data and scaled to the 5.1 cm airfoil. 

The frequency was scaled by plotting the autospectra of the 

longitudinal component, of turbulence for each case on a common 
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scale, then overlaying the spectra. Alignment which gave the 
best agreement corresponded to a frequency scaling factor of 10. 
Since the scaling was being done from the fully turbulent region 
of the large airfoil to the transitional region of the small 
one, this overlaying scaling method was felt to be the most 
valid one. The actual values of the correlation lengths were 
taken from Figure 46 where the Reynolds number effect was in- 
cluded. The scaled lengths were then used in the prediction. 

Results of the predicted spectrtim for sound radiation nor- 
mal to the airfoil chord are shown in comparison with the mea- 
sured data in Figure 46. Agreement is good considering the 
fact that the correlation length dependence on frequency was 
assumed to be scalable to the transitional flow regi.me to the 
fully turbulent regime. 

It is important to observe that the roll off in the rad- 
iated SPL at higher frequency is predicted by the model. This 
effect is critical in predicting useful spectra since any sub- 
jective evaluation will be very sensitive to the high end of 
the spectrum. The importance of this can be seen in the sam- 
ple rotor case to be discussed later in this section. 

A radiation pattern for the airfoil is shown in Figure 47. 
Radiation patterns for different frequencies are plotted in 
the figure. A substantial difference can be seen to exist with 
the levels for different frequencies. The differences can be 
interpreted by envisioning two dipoles, one aligned parallel 
to the airfoil chord and the other perpendicular to it. The 
parallel dipole has magnitude corresponding to the longitudinal 
component of turbulence and the perpendicular dipole corresponds 
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to the transverse component. In those parts of the spectrum 
where the power spectral level of the transverse component is 
greater than the longitudinal component, the dipole in the 
perpendicular direction will dominate. The observed dominance 
of the transverse component of turbulence in the transitional 
range of Reynolds number may explain the apparent single dipole 
nature of the acoustic data taken there. As an example, the 
data shown in Figure 48 illustrate this tendency. Notice even 

here there is a greater fullness to the measured pattern than 

in the theoretical dipole. This fullness may be explainable 
in terms of the reduced, but still existant, contribution of 
the longitudinal component at this Reynolds number. Turbulence 
data taken at the time indicated the normal component of tur- 
bulence to be about 10 times the level of the longitudinal com- 

ponent (18). 

When the cross terms are included the radiation pattern 
will no longer be symmetric but will be distorted in the forward 
and aft arc depending upon the net shear stress in the wake. 

A net shear stress will result when the wake is non- symmetric ; 
this would occur for highly loaded airfoils. 

6.2 Scaling 

Before discussing any detailed scaling method it is impor- 
tant to be aware of the existence of the two flow regimes which 
may be present in any particular case. It appears that laminar 
separation induced vortex shedding will contribute in a major 
way to sound generation at transitional Reynolds numbers. This 
mechanism is characterized by pure tones or nearly pure tones 
corresponding to the frequencies in the vortex shedding. 














103 


Patterson et al (17) have made extensive measurements of this 
effect and were able, through a mapping over a wide range of 
Reynolds number and incidence angle, to estimate the flow re- 
gimes where the tone mechanism would be present, A typical 
plot from their findings for an NACA 0012 airfoil is shown in 
Figure 49. The tone appears to be present at higher Reynolds 
numbers when the airfoil is operated at non-zero incidence an- 
gles. This author has observed similar effects (16); not with 
the pure tones observed above but rather a narrow band random 
type of spectrum. As the Re 5 molds number is increased the tone 
tends to change into a narrow band random shape and eventually 
vanish into the wide band turbulence spectrum, 

A scaling procedure was suggested for frequency in refer- 
ence 17 which used a relationship for the Strouhal number based 
upon boundary layer thickness. Recall the form 
St - 2 f U 

and the laminar boundary layer thickness was written 


Slam = K C 


-.5 


This resulted in a frequency scaling dependence of the form 

f oc / (C'^^i *^) 

The scaling in this way was used successfully to scale the tones 
they observed in their experiments and also the peaks reported 
in the spectra in reference 16. 

An approach similar to that outlined above was used in an 
attempt to obtain a frequency scaling rule for the fully turbu- 
lent region. 

In what follows a limitation to flows where the tone does 


not exist is made. This generally will mean chord Reynolds 


Incidence Angle, deg 


® Vortex Shedding Noise Observed 
^ No Measurable Vortex Shedding Noise 


Note: This figure taken 
from ref . 1 7 


Chord Reynolds Number 


Figure 49 Observed Regions of Vortex Shedding for a 
NACA 0012 Airfoil as Function of Incidence 
Angle and Reynolds Number 
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greater than 3 x 10 . This restriction will have little prac" 
tical significance since most machines of interest operate at 
the higher Reynolds numbers. 

The scaling procedure will be established for the corre- 
lation lengths, the frequency, and the mean flow parameters. 

This present method is suggested for use only with small values 
of Incidence angle since significant amounts of non-sero inci- 
dence angle data are not yet available. The method would be 
expected to underestimate the actual spectral levels for the 
more highly loaded airfoils. 

6.2.1 Correlation Lengths 

Correlation lengths are ssumed to be proportional to the 
turbulent boundary layer thickness. Dependence of the boundary 
layer thickness on Reynolds number for fully developed turbulent 
flow on a flat plate has been shown to be 
S^/ C ocRe"*^ 

This approximate dependence was also observed for the correla- 
tion lengths measured for the 25.4 and 91.4 cm airfoils. Based 
upon this, the scaling can be written 

(l^/C) = (i^/C)g]^ <Re/Regj^)'-^ 

where the Reynolds, number for the large airfoil is Regj^=2xl0 . 
The above result is used as a frequency independent scaling 
parameter which means it is used across the entire spectrum with 
no further modification. 

6.2.2 Frequency 

An approach similar to the one described for the tone 
'frequency scaling will be outlined here for the turbulent flow 
regime. Defining the Strouhal number in terms of the turbulent 
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boundary layer thickness 
= 2 f / U 

and assuming there exists a common Strouhal number for the en- 
tire turbulent regime, this relationship can be used to scale 
the frequency. Dependence of the boundary layer thickness on 
the chord Reynolds number can be expressed in the same way as 
was done for the correlation lengths. Combining with the 
Strouhal number yields a scaling factor for frequency which is 
f = fgj (U / (C / 

with: Ug^ =45.7 m/sec and = 91.4 cm. 

One opportunity exists to check the scaling rules for 
length and frequency by using the measured data for the 25.4 
cm chord airfoil. Since only data were taken for the longitu- 
dinal component of turbulence, comparisons are limited to this 
component. The spectra of the correlation lengths for the 91.4 
cm airfoil were scaled to the 25.4 cm case using the above rules. 
The scaled lengths could then be compared with those actually 
measured. Pertinent variables for the scaling are: 

91.4 cm airfoil 25 .4 cm airfoil 

U = 45,7 m/sec U = 23.7 m/sec 

Re = 2 X 10^ Re = 2 X 10^ 

incidence = 0”^ incidence = 0° 

Results are shown in comparison with experimental lengths 
in Figure 42. The frequency scaling appears to be fairly good 
since the shapes of the scaled and measured spectra are roughly 
the same. The actual values of the lengths are scaling with 
a maximum difference of about 20 to 30%. However, it is appar- 
ent that the scaling is different in different parts of the 
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spectrum since the shapes of the curves are not the same. As 
a first approximation the procedure doesn’t look too bad. 

6.2.3 Mean Turbulence Levels 

It is recommended that the actual results for the 91 cm 
airfoil as non- dimens ionalized by the free stream velocity be 
used directly with no further correction. Measured RMS turbu- 
lence levels were very nearly the same for all three airfoils 
so the usage of the data from the large airfoil should introduce 
little error. 

6.2.4 A simplified Sound Radiation Model 

A simplified form of the prediction model may be establish- 
ed for use in trade studies where detailed aerodynamic data are 
not available. 

Ignoring the effects of the cross terms and the spatial 
distribution of the wake properties, the model can be expressed 
in terms of mean wake parameters and a simple algebraic form 
written. The power spectrum for the radiated sound density: 

UfU./ (iS) 

•f- Ij ^3 , ) 


where it can be seen that the dependence on the velocity is to 
the sixth power as expected from the dipole force mechanism. 
This sixth power results from the second power directly on the 
velocity and the fourth from the frequency which is itself 
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approximately proportional to the velocity. The sound is pro- 
portional to the wake thickness the airfoil span S, the 
stream correlation length L^, and the correlation volume approx- 
imated by the product Lj^L 2 L^. The velocity used probably should 
be a mean velocity in the wake; a value of 0.8 U would be rea- 
sonable in view of the measured velocity obtained in the exper- 
imental data. 

The above method will be expected to work best for air- 
foils with small lift coefficients. For more highly loaded air- 
foils which better characterize real machines, one can expect 
the prediction to be low. It is expected that incidence angle 
correlations will be possible once more data are available, 
particularly on the effect of the incidence angle on the span 
correlation length. 

6.3 Rotor Results 

The scaling method established in section 6.2 was used to 
supply information to enable the prediction of the spectrum of 
sound generation by a subsonic rotor. The design chosen was 
to roughly approximate what might be expected for a large sub- 
sonic fan as may be used in an STOL type of aircraft. The de- 
sign was approximated by a single source for each blade in the 
rotor. This is a crude approximation but in view of the amount 
of computer time required to carry out a distributed source re- 
presentation and the absence of any particular design it was 
felt sufficient to demonstrate the effects with a simple single 
force per blade. The rotor program has been devised to predict 
the sound for the direct problem of a given aerodynamic design 
and must be used in this way. 
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The design details for the particular rotor to be studied 

are as follows: 

Diameter; 4.56 m (15 ft) 

Chord: .91 m (3 ft) 

t/C: 0.09 

Tip speed: 182 m/sec (600 ft/sec) 

Source Radius: 0.8 

Stagger angle: 65° 

Blade number: 12 

Using these variables and the scaling procedure, the predicted 
sound radiation spectrum for three radiation angles measured 
from the inlet rotor axis are shown in Figure 50, There are 
some differences in the pattern at the different angles but the 
overall effect appears to be to radiate essentially the same 
amount in all directions. This essentially uniform radiation 
pattern has been observed experimentally by the author for 
several rotors ranging from 6 inch diameter compressor rigs 
all the way to the Pratt and Whitney JT9D engine when operated 
in the subsonic Mach regime. 

Effects of source motion on the radiated sound can be ob- 
served in a simple way using the rotor model. If the RPM for 
the rotor is set to zero but the relative velocity onto the 
blades is left at the same value (both are Independent inputs 
to the computer program) the rotor behaves as a stator. This 
stator configuration can then be used to illustrate the effects 
on the radiated sound spectra by Doppler changes. The two spec- 
tra for the rotor and the stator are shown in Figure 51. There 
is very little change in the radiated sound spectrum due to the 
motion of the sources. This is to be expected since the shift 


Radiated SPL at 152 


Angle to Rotor Inlet 

O 0“ □ 45° 


O 90° 


Note; Levels are 6% bandwidth 








F r eqnency , kH z 


PNLgQO = 103,5 PNdB 


Figure 50 Radiated SPL - Spectra for Wake Related Sound by Large 
Subsonic Rotor - Dia. = 4.6 m, = 183 m/sec 



Frequency , kHz 


Figure 51 Illustration of the Effect of Rotation on the Radiated SPL, 
Design are Identical, One Source Set is Rotating 
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in the spectnim is to neighboring levels and since the spectra 
of the two point velocity correlations are essentially flat 
there would be little effect. However, if the source spectra 
were more narrow band in nature as for the case where vortex 
shedding might be present, then the effects of the source mo- 
tion could be substantial. 

The last point to be discussed is the subjective sound 
level of this rotor. The 500 ft ( 152 m) sideline noise level 
as indicated by the Perceived Noise Scale was computed to be 
103.5 PNdB. This is the noise due to the wake related broad 
band effects. Any contribution by the pure tones or due to 
inflow turbulence would add to this figure. This is a fairly 
high level. It results because the frequency of the sound hap- 
pens to fall at levels where the htiraan ear is most sensitive. 
There are currently proposed engine designs which are similar 
to the one used here and they may have this noise problem. The 
idea of these designs is to use a relatively large number of 
rotor blades since cancelation effects tend to reduce pure tone 
noise. Based upon the present results, the trade to wake noise 
may be too much. The quality of the noise will be much like 
the jet noise from current turbojet engines but with the fre- 
quency shifted to 4 or 5 kH 2 from the present 1 kHz. The sound 
would be a scream'. 

These large fans may be amenable to design procedures which 
could reduce the wake related noise. If boundary layer control 
could be used on the blading, the reduction of the correlation 
volumes in the wake mighh result in reduced noise generation. 

A great deal of research would be needed to develop a viable 
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control system for a rotor of this size and complexity, 
would envision a very long and costly development period 


One 
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7.0 Conclusions 

1. It appears based upon the agreement obtained between mea- 
sured and predicted sound spectra that a sound prediction model 
based on approximating the force fluctuation on an airfoil by 
the fluctuation of momentum in the near wake is a useful method. 
However, to build confidence in this approach more experimental 
sound radiation data from single airfoils is needed. 

2. The effect of the Re 5 molds number can cause large variations 
in the extent of correlation in the span direction in the wake. 
Between the transitional and fully turbulent flow regimes 
changes in the span correlation length of an order of magnitude 
were observed. 

3. A major effect of the incidence angle appears to be in the 
normal correlation length. A change of 6° resulted in a factor 

of nearly 4 increase in the correlation length in this direction 
with significantly less in the other directions. 

4. Using the results for the single airfoil model a scaling 
procedure was developed and used to predict correlation lengths 
at another Reynolds number with fair success. More work is 
needed before an accurate scaling method can be expected. 

5. When extended to a rotor configuration it was found that 
large scale subsonic fan rotors may be expected to generate 
wake related noise of significant subjective levels. The possi- 
bility exists for reducing noise by these fans through boundary 
layer control on the rotor blading. Development of this kind 

of flow control would be expected to require a substantial 
investment of time and money. 
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Model for the Generation of Aerodynamic Sound 
by a Single Airfoil 

The model proposed for aerodynamic sound generation by a 
single airfoil is based on Curie's (5) solution of the non- 
homogeneous wave equation first described by Lighthill (2), 
Curie's solution is reduced for the case of an airfoil opera- 
ting in a smooth flow. The only turbulence present in the flow 
is that in the wake of the airfoil. The solution is reduced 
to a form which relates the spectrum of sound generation to the 
power spectra of the two point velocity correlations in the 
wake. An outline of the wave equation, its solution, and reduc- 
tion is presented in this appendix. 

Non-homogeneous Wave Equation and General Solution 

This wave equation is formed by combining the equations 
of continuity of mass with those for conservation of momentum. 
The continuity equation used is 

and the momentum equation is used in the form 




dXj “V 


For the energy equation, the wave propagation is assumed to 
be i sen tropic so that the density and pressure may be related 


by 


^ (-p %) - % (p Stj) 


where a^ is an average wave speed in the radiative medium (air) 
The procedure is to take the partial derivative of the continu- 
ity equation with time and the partial space derivative of the 
momentum equation and subtracting the two results, A form of 
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the wave equation may be written 

2 -s2 2 


I# ~ '^dXidx^y ' 

If the isentropic relationship is used and the viscous stresses 
are neglected as small compared with the turbulent stresses, 
the wave equation takes the simpler form 


Cou-u) 


For an unbounded medium the solution is given by the Kirchhoff 
retarded potential solution (Sobelev, 4) but a more general 
form is needed for the case where bodies are immersed in the 
flow. A form of this type is also given by Sobelev but it is 
the reduction of the general solution to a form pertinent to 
sound generation which was the contribution of Curie (5). 

Curie has developed a solution which may be written for station- 
ary surfaces as 



The delta function reflects the retarded time at which the in- 
tegrands must be evaluated. The retarded time is the time it 
takes a disturbance to propagate from the source to the receiver. 
The integral must be evaluated at this retarded time unless it 
can be shown that negligible error results with its neglect. 
Fortunately for many cases in turbulence related noise generation 


are 
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the retarded time effect may be safely neglected. These 
typically low Mach number flows. 

Before considering the retarded time some other reductions 
are possible. Details of these reductions have been presented 
by Curie however for convenience and continuity they will be 
outlined here. 

Of the two terms in the solution the stress term is the 
more complicated so consider it here. Define r = [ x - y\ 
and write 







Now repeat the same procedure with the second derivative and 




becomes 


f f dr^_rn d%(y.t‘) 2 dTa(tf) 2 -r t')l 

S(t'-t -i- -a J 


Since 



it can be shown that 


t* , ' r 


In order to get a measure of the relative magnitudes of the 
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three terms in the integral the method outlined in ref. 19 
is used. 

For data of some typical frequency f we can write an equi- 
valent wavelength from the expression 

•f X = 0*0 

With this the three terras have the relative magnitudes 

±n£ -r- 2 ^ 

and in terms of the wavelength 

So that for ()s/r)-i<l the second two terms may be neglected and 
the solution written 



In an exactly parallel way, the surface term may be reduced 
to the form 

A further reduction may now be made if the origin of the coor- 
dinate system is located within the active turbulent region 
(the airfoil wake). Then the source dimensions are small com- 
pared with the observer coordinates, y^« x^. 
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The far field solution may now be written 



This forms the far field solution to the wave equation with 
immersed bodies. Originally it was planned to include the 
volume stress terms in the solution and attempt to measure 
these as well as the force terms. This was abandoned for the 
reasons given in the body of the report that the dipole terms 
appear to dominate this mechanism of sound generation and there- 
fore the volume terms would be of smaller importance in the 
sound generation problem being considered. If considerations 
were being made of higher Mach number flows there is evidence 
that the "quadrapole" terras may be of significant importance. 
Since the present experiment was limited to low Mach numbers, 
there appears to be little justification for the inclusion of 
those terms in this analysis. For these reasons, the rest of 
this development will be directed to the surface terms. It 
should be pointed out that a great effort has been extended in 
developing an understanding of free turbulent flows since these 
are the source of the very critical jet noise problem facing 
the commercial airplane industry at the present time. The pro- 
blem is very complex but substantial progress has been made and 
it is expected that these same techniques will be applicable 
to the wake turbulence mechanism. 
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Also, it should be pointed out that sufficient data have been 
obtained in the present study to allow the establishment of 
the forcing function for the voliime tern. A very large data 
reduction effort would be required to reduce the data to a 
form useable for the volume source specification. These data 
are on magnetic tape and such a reduction effort may be under- 
taken in the future. For the present purposes it is felt that 
efforts may be more profitably directed to the study of the 
dominant mechanism of the fluctuating force on the airfoil. 
That alone has proven perplexing enough. 

Conditions for Neglect of the Retarded Time 

The clearest presentation, at least for this author, of 
the implications of the neglect of the retarded time in sto- 
chastic processes has been presented by Morse and Ingard (10). 
It is this general idea which is paralleled below for the pre- 
sent case. For simplicity in the development of the idea con- 
sider a function of form similar to the solution of the wave 
equation. Write 





Form an autocorrelation using the expectation as 


=E) 


/ r 




(I Ky,t')6(t'-t+^) P(?,t%x)s(H:+t-b+^cl2dyclt'| 

E [p(^.t') P(€,t'+t)j S(l:-t S (t'+i:- 1 didfdt' 
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The process is assumed statistically stationary so that the 
only significance is the time difference between the two signals 
and not the absolute time. This difference between the times 
of the two signals may be written 


(t+xV-t = t+X- 


= X + 


\t-vl lY-tl 


do 

Expand these functions to 


do 


'It 


\/(x.-y.)2 = (lXl^ + lyl*-2x.y,-2)f4ya-2Jf3/5} 


- A - 




/Xl 


1x1' 


m 


To the first order in y/x the expansion may be written 

i$-4\ = 

J. 

w 

and similarly for the other source position vector, z: 


\i-t\~ \i\- 


y * fe 


Using these results the time difference becomes 

tto a-o a.olx\ 

Notice that *x/\x\ is the unit vector in the direction of the 
observer. The second group is the time for a sound signal to 
travel from "y to "z in the direction of the observer. 

Write the correlation in these terms as 


1 
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f 


> 


= E -TpCv) P(t,x+-Sx) • dldy 

■'w ^ 


It can be seen that if the additional time caused by results 
in no significant difference in the correlation between the two 
functions P then we may neglect the retarded time difference S'z 
Consider this further. Notice: 


- y) = 1^ = a projected length of correlation 
in the observer direction 


Therefore, write: 

X + St: - "C + 

and consider it In terms of the Fourier transform. The argument 
for the transform will be 

co(x + B-c) = -f- = loz -f- ?i0£. 

It can be seen from this that there will occur phase changes 
which are dependent on the ratio of the correlation length to 
:.he wavelength of the sound. Some interference will occur while 
the sound traverses the correlation length if the wavelength 
is not much larger than the correlation length. A condition 
for the neglect of retarded times can be written as 

For the single airfoil data this requirement appears to be 
satisfied. However for the higher Mach numbers it may not be 
sufficient over the entire frequency range of interest; this 
must be remembered when extending these results to higher Mach 
number flows. If the retarded time is neglected the solution 
may be written in the form 


£.ir Jtc 
A 
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ay 

This is the form which will be reduced using the momentum equa- 
tion so that the surface integral can be replaced by a volume 
integral over the wake near the airfoil trailing edge. 

Change of Surface Pressure Term to Volume Momenttim Term 

Because of the difficulties in making static pressure 
correlation measurements on the surface of airfoils, the pre- 
sent approach of relating the force fluctuations on the surface 
jf an airfoil to the momentum fluctuations in the wake was 
taken. A more detailed discussion of how this approach was de- 
cided upon and its limitations are presented in the body of the 
report. Details of the development of the actual model are 
presented in this section. 

What will be done is to relate the surface static pressure 
on the airfoil to the gradient of the static pressure in the 
wake volume using the Gauss divergence theorem. The resulting 
(grad p) term will be expressed in terms of velocities in the 
wake through the linear momentum balance equation. Covariances 
of the resulting terms will provide the relationship between 
sound generation and the spectra of the two point velocity 
correlations in the airfoil wake. 

First, consider the surface integral over the boundary of 
the control volume shown in figure 55. The control volume is 
such that the only non- steady velocities (other than lower 
order acoustic velocities) are confined to the turbulence of 


n 
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Figure 52 Control Volume Used in Application of Momentum Equation 
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the airfoil wake. For the simply connected region shown the 
divergence theorem can be used to write 



Taking the derivative of the expression allows the regions 
having steady properties to be eliminated from the integrations 
and simplifies the algebra. The result is 



AUi^dSCf) + 

aiHoil 


UIBKS BXIT PUIME 


Through the divergence theorem these are then equal to the 
unsteady voltime integral of the (grad p) in the wake. Therefore 
write 



The (grad p) integral may now be related to the velocity fluc- 
tuations by means of the linear momentum equation. For our 
case of incompressible flow 

-- grad 

so that 


Jy 

The second term may be expressed as a surface integral by again 
using the divergence theorem. The total expression becomes 




P 
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The force on the body is to be approx? niated by the momentum 
fluctuation in the wake immediately downstream of the airfoil 
trailing edge. The assumption is that the force on the blade 
can be approximated by a kind of correlation volume of momen- 
tum nearest the trailing edge. This assumption then eliminates 
the dependence upon the terms in the exit plane integral since 
these will be uncorrelated and may be ignored. 

Using these conditions the relationship between surface 
pressure fluctuations and those of wake momentum become 

ourh! ^ 

with this equality the radiated sound can be written. It is 
necessary to first introduce a negative sign since the sound is 
generated by the force on the fluid where in the above the force 
of the^luid on the body has been determined. Write for the 
radiated density 

, . ■"/?> Ji ^ f ^ 

pU,t)'A = itl^ dt ^ 

Notice that the density has been taken outside the integral 
under the assumption of incompressible flow. At this point an 
examination of the consequences of this assumption and the 
expected error due to density fluctuations will be made. 

Estimate of Error Due to Neglect of Density Fluctuations 

As a means of getting a measure of the error due to the 
neglect of the density fluctuations the energy equation for 
one dimensional compressible flow will be used. 
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Write the energy equation 

and introduce the normal turbulence assumptions 

UL-U-t-u! 

Factor the energy equation into the form 

ZCp9-(^ “ ‘^) 

Where the fluctuations in static pressure have been assumed 
small compared with those of velocity and density. This can 
be taken as generally true for free flows such as the wake- 
Rewrite the energy equation as 

/oo^ = Zep9i;(/o-/o^) 

and insert the turbulent velocities and densities to yield 

(p+/3')(0 +ZUtl■^■u.'u') = ZCf=dt(/° 

Solve this for the density fluctuation: 

XUp ^ ~ ^ — ZCp&^^ 

For air (^= 1.4) and the fluctuations in density are related 
to those for the longitudinal component of turbulence by 



This suggests that for M < 1 the density fluctuations may be 
safely neglected. This is therefore within the context of the 
other assumptions made in the current model. 


131 

Sound Radiation in Terms of Two Point Velocity Correlations 
With the assumption of constant density the covariance 
of the radiated sound density may be simply related to the co- 
variances of the two point velocity correlations. 

Separate the velocity into a steady part and a fluctuating 
part in the typical turbulent assumption: 

u-i ( Y,t) = Ui ty) + (%t) 

and use the direction cosines 

A 

Refer to figure 4 and write the direction cosines to agree 
with the nomenclature of the figure. They become 

-a. ^ 

^ O vrt X * ^7 VS) 

—3— = = COS T 

Further, assume the velocities arc represented by the following 
components ; 


Using index, notation, write the solution as 






^ 7 ra|/x/ 


A 

n 






then expand and normalize with the free stream velocity to gel 


peuypo^ 


PqU. 

4-jtaliii 
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Define the following coefficients to help simplify the formu- 
lation of the covariances: 


A t^.-^,eosY 

^Z. 4-Tra^/xl 


Write the expected value of the radiated sound density as 
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Replace the expectations with the cross covarianc.' functions 
to obtain the form 

•'y I I « t. 




f 
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In order to cast the solution into a form usable with the exper- 
imental data the Fourier transform of the covariances must be 
taken. The one-sided power spectral form of the transfrom will 

be used. This may be written 

. + 00 

X ^ -s. X -Loot / 

Z R,y(x,z)e dz 

Loo 


Using this and the identity 




uu 


write the solution in terms of the cross power spectral densi- 
ties of the two point velocity correlations. The result is 




^ \ 
y,0) 


This is the general form of the present solution for sound 
radiation from a single airfoil in a smooth flow. It is the 
basis for the model and these are the quantities which were 
measured in the experimental program. 

Determination of the Momentum Fluctuation 

For the rest of this appendix the summation convention is 
suspended. The spectrum of the radiated sound can be written 
in the index form as 


z z 



L = f 
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where 



Now change variables and let one of the position vectors of 

the sonrces be defined relative to the other as 

■ ^ 

2 = y +>1. 


and write 






r/. Y2> oV/?/ ^y, 


where: u^u . (y, ,y2,y3 ,dJ ) = the distribution of the auto- 

^ spectrum in the wake 

) = a coherency spectrum 

The term labeled a coherency spectrum in the above is similar 
to the normal coherency as used in texts on time series analy- 
sis with the difference that the normalization is done with 
only one autospectrum, that of the signal on the stationary 
probe, and not with both signals as is the normal procedure. 
With this type of a definition it is possible to have values of 
the coherency greater than one. The Schwartz inequality does 
not apply. 

Notice that the coherency spectrum has been assumed inde- 
pendent of y2 and y^ but the dependency upon y^^ has been re- 
tained. The argument to support this is that the magnitude in 


the 72 direction decreases so rapidly (as one nears the edge 
of the wake) that differences in the spatial distribution of 
velocity will be of small effect in the overall Integrated re- 
sult. In the y^ (span) direction the assumption is valid since 
spatial homogeneity may be correctly assumed in this direction. 
The fact that the stream dependence cannot be negJected is due 
to the fact that significant differences may exist in this di- 
rection. 

Usage of Correlation Length Concept 

It will prove convenient in terms of handling the experi- 
mental results to carry out the integration of the coherency 
function in each of the coordinate directions at each frequency. 
The lengths that result will be termed correlation lengths as 
a method of labeling. Notice that implicit in this method is 
the assumption that the integration can be done in each direc- 
tion independently of the others and the integral can be ex- 
pressed as a product of the three lengths. 

Each length becomes a function of and 0 > . The integral, 
may bo written 

i " 


where the correlation lengths may be different for the various 
two point velocity correlations; hence the dependence on 
is retained In the lengths. 

The sound radiation due to one particular u^u^ may be 


written as 
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C^i^j 3 ^ Yi j ^") L^ O^i^jj Y}y**0 C^^jj Yi > 

y» 

“j/ ^i^j ^■(%'f%>Si,'^)<^V3^yi^y, 

% % 

In line with the assumptions of the model, assume the relevant 
volume of integration to be that correlation volume immediately 
do\mstream of the airfoil, is then taken over the range -L^ 

to and the expression for the sound radiation is 

L-=i /^/ 

<ii^(yi,y3,‘^)dy3°fYz 

This is the final result for the model which will be used to 
relate sound generation to wake turbulence. The first integral 
(that of the longitudinal component of turbulence) will corre- 
spond to a dipole oriented along the chord direction. The la=t 
integral, the transverse component, will be a dipole oriented 
along a line perpendicular to the chord of the airfoil. This is 
the dipole normally used in models for the prediction of vortex 
shedding noise. The cross terms will act to modify the radia- 
tion pattern from the symmetric form. If, however, the wake is 
symmetric these terms will integrate to zero. 
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APPENDIX B 
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Wake Related Sound Generation by a Rotor 

The mathematical model used to predict sound generation 
by a rotor using the general results of the single airfoil 
was that of Ffowcs Williams and Hawkings(21) . There is avail- 
able another analysis based on the Lowson theory for moving 
sources which was recently presented by Morfey and Tanna(22), 
The quantity of algebraic detail required to extend this latter 
analysis to include both forces and stresses is considerable 
and the application to specific designs was felt to be sub- 
stantially more laborious than the method of ref, 21. Actually 
there are many rotor analyses springing up like mushrooms but 
these do not appear to offer any particular advantages over 
those already mentioned for practical designs. 

An outline of the original analysis by Ffowcs Williams 
and Hawkings will be given in this appendix with the bulk of 
the space used to develop their results into a form usable for 
the particular problem of wake related sound generation. For 
a general discussion of multipole analysis, of which this is 
an example, the reader is referred to Doak(6) and to Morse and 
Ingard( 20) . 

The sound generation is taken as governed by a non- 
homogeneous wave equation of the form 

- al V/5 = <3 (x,t) 

vrhich has the solution 

s(t'- 1 >y&o)cly 
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and in terms of its generalized Fourier transform, 

As may be seen from Curie’s solution to the Lighthill equation, 
the forcing function may contain simple sources, dipoles and 
quadrapole elements. In general the form is 

- s(y,t) ^ ^ 

3yi ayj 

So that Sj D. , and T- • are the source strength densities of 

X ^1 

the three most common source types. For present purposes, only 
the force and stress densities will be considered. 

The source expression may be Fourier transformed to ex- 
press it in the form 

= f f Q(f,t)<s. didf 

~'yJt 

from which, 

Q(k,c^ = S(k:,cS) i- L kiO^Ck,(S) - 


The components radiating in the direction of the observer can 

ri _ CjJ ^ 

be taken as the spectral level with, K — ~ 


thU: 




The general solution was then developed by substitution and 
reduction of the above equation inserted into the general 
solution as it was originally put down. What resulted were 
the following forms for sources moving in a plane circular 
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motion (a rotor). Details of this reduction are lengthy and 
best left to the original paper. The results for the two source 
types of interest here, i.e. dipoles and quadrapoles, are pre- 
sented below. 


ROTATING FORCES (DIPOLES) 

-in(ir/2-eie) 
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ROTATING STRESSES (QUADRAPOLES) 
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These solutions are written for a point dipole having a 
generalized spectrum d^(0f^,«^) and a point.' quadrapole with a 
generalized spectrum t^j (0j^,o(.) . The spectrum is the Fourier 
transform of the time dependent source strength. The position 
is the angular location (8]^) of the kth source at time, t - 0. 
XL is the angular frequency of the source and <f> is the radia- 
tion angle measured from the inlet fan axis. 

In order to put these solutions in a form usable for the 
prediction of wake related sound, an auto power spectrum of 
the radiated sound density must be formed. This will result 
in cross power spectra of the various terms in the solution. 
The details of this procedure for combined forces and stresses 
are developed in the following paragraphs. 

Reduction of the Solution for Multlbladed Rotors 

The total radiated sound is taken as the sum of the two 
components written above. 


For one blade of a rotor, write 


The total for the entire rotor is the sum over the number of 
rotor blades. 

B B 

/oat, CO) = ^ ^ 

fe=/ 


To obtain the radiated sound power spectrum multiply the above 

equation by its complex conjugate to form the modulus squared. 

B 

\p{t,io)f= [E /°d,le 


te-i 








expand the form to 
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l/0(X,Cb)l‘ 



-}f 


i-\ 

For algebraic convenience consider the first term in the 
following development as representative of the others. 

g ^ ^ ■Hco 

fiL 

' o m=-o= 

S . 44- 


fcr fei" 

-Lnfir/g-fflfe) -i-mOr/Z-Sji) f eoR^sm^\-r ( coR^sm<A 

e e j„[- /''"'i ao J 

where d^ is the component in the direction of the observer. 
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Consequences of the Assumption of Stationarity 

If the force fluctuations are assumed to be statistically 
stationary, then a reduction in the complexities will result. 
For stationary random processes the expected value of a form 
similar to the solution may be expressed. 




where S^y.(£) is the cross power spectrum. This states that for 
stationary processes the Fourier coefficients are uncorrelated. 

Using the above result the only terms which will remain 
in the summation are those for which m = n. The solution re- 
duces for this case to 



The Assumption of Uncorrelated Blades 

If the process at each blade is assumed uncorrelated with 
those at each of the other blades, then the double sum on the 
blade number may be reduced to a single sum. 

If the processes are mutually uncorrelated from blade to 


blade 






rf 


Notice if a model is being developed for inflow turbulence 
related noise generation the reduction made here would not be 
correct since the scale of the turbulence would be expected 
to encompass more than one blade in the interaction and 


response. 
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For the present case of uncorrelated blades it is sufficient 
to write the dipole solution. 

j 2. B 

CO 

f CO Rs Sw4>\ 

V ■““a; J 


where the identity has been used. 

Consideration of the other three terms in the solution 
will exactly parallel the above arguements so it is sufficient 
just to write the result. The solution becomes 

IpCXjOi^l " ZZ Z— 1 

V\=-co te.= i ^ 



Notice that there appears to be an imaginary part to the 
modulus. This of course is not possible. Actually the appar- 
ent complex nature of the expression is not really there since 
the non-vanishing part of the term in square brackets is the 
imaginary part which is multiplied by - T and produces a real 
number. This may be shown by expressing the result in terms 
of co“ and quad-spectral components; the co-spectral components 
will cancel leaving only the quad components (see Bendat and 
Piersol, page 83). 
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Representation in Terms of Wake Parameters 

Since the rotor model is developed here for point sources, 
the distributed results to be taken from the single airfoil 
results must be expressed as an equivalent point dipole. The 
method used for the force will be developed. The stress term 
was not developed farther for the wake case since for the range 
of variables studied here the sound generation was dominated 
by the dipole mechanism. Originally it was expected that 
the quadrapole terms in the wake would be included so the rotor 
analysis was set up to handle both. It is still expected that 
for higher Mach numbers wake stress terms will assume Importance. 
The rest of this rotor development will consider only the force 
fluctuations on the blading. 

The total force on the blade is taken as the integral of 
the static pressure over the surface of the blade. 

5 

As was seen for the case of the single airfoil, this static 
pressure distribution could be written in terms of the momentum 
fluctuations in the near wake of the airfoil. 

dSC'J') = j |:^(.paOciy 

In terms of frequency the force may be expressed 

cUe^.io) = \PTto^ j Ui,^Y,w)dy 

y 

This is the form used in the rotor model. 

Substitution of the force into the general solution for 
the sound radiation gives the force modulus for the component 
in the direction of the observer. 








write the force modulus 
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The computer program which was developed to predict the rotor 
sound included both the stress terms and the force terms as 
v^ell as the cross terms for the correlation between the stress 
and force. As the program is used for this current work only 
the force, or dipole, part will be employed. It is expected 
that as more understanding of the relationship of the contri- 
bution of the wake stress is developed these results can be 
included in the rotor computations. 


Nomenclature 


dj.(6k,f) 


E 

f 


tT'V 


u,u . 
. 1 J 


M 

n 


P 

Pk 

Q(x, t) 


A 1 J 
’^k 

S(y, t) 

s 



Tij(y.t) 


speed of sound in air (m/sec) 

specific heat at constant pressure (Kcal/Kg °K) 

dipole rotor source term (N/m ) 

spectral level of blade force component 
in the direction of the observer (N) 

base vector, rotor model 

expectation operation (expected value) 

frequency (Hz) 

2 6 

density power spectrum at the observer (Kg /m^Hz) 
base vector, rotor model 

one sided power spectral density (m /sec Hz) 
base vector, rotor model 
typical correlation length (m) 

correlation length for the k coordinate direction (cm) 
Mach number 

unit normal to control volume 

2 

static pressure ''N/m ) 

2 

force per unit a”ea (N/m ) 

total rotor source term (N/m ) 

^x - y\ distance from source to observer (m) 

2 2 

cross covariance function (m /sec ) 
direction cosine, X|^/|x\ 
simple source term, rotor model (Kg/m^) 
airfoil span (m) 

2 

Lighthill stress tensor (N/m ) 

(used in single airfoil model) 

quadrapole rotor source term (N/m ) 
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,f) spectral level of the stress coinpo-gent in 
the direction of the observer (N/m^) 

time (sec) 

free stream velocity (m/sec) 
mean velocity (m/sec) 

fluctuating component of velocity (m/sec) 
two point velocity correlation (ra /sec ) 
position vector for sound observer (m) 
position vector for source (m) 
position vector for source (m) 
density of fluid (Kg/m^) 
radian frequency (rad/sec) 
wake thickness (cm) 

angle from perpendicular to airfoil chord 
to the sound observer (deg) 

2 

shear stress tensor (N/m ) 
wavelength (l/m) 

arbitrary function, used in Appendix A 
ratio of specific heats ( JT = 1.4) 
total temperature (°K) 
static temperature (°K) 
fluctuating density (Kg/m ) 
displacement vector, "z - "y 
rotor frequency (1/sec) 

radiation angle from inlet rotor axis to 
sound observer (deg) 

angular location of rotor source at t = 0 


(deg) 


